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ABSTRACT 


\ 

This  dieeetiteafc4ett'  presents  the  results  obtained  in  an  experimental 
and  theoretical  study  of  the  scattering  of  underwater  acoustic  waves  from 
a rough  surface .'■j  Four  rough  surfaces  J2  in.  by  32  in.  in  size  were  made 
of  a pressure-release  material.  The  reliefs  of  these  surfaces  were  approxi- 
mately Gaussian /and  vere  characterized  by  an  exponential  covariance  function. 
The  amplitude  ajid  relative  phase  of  a pulsed  wave  scattered  in  the  specular 
direction  were  pleasured  at  a frequency  of  95 3 kHz  for  grazing  angles  from 
6 deg  to  80  de^.  Both  the  directional  source  with  a beamwidth  of  9 deg 
measured  at  thel  half-power  points  and  the  omnidirectional  receiver  were 
placed  at  a slant-line  distance  of  60  in.  from  the  scattering  surface  for 
all  measurements . measurements  were  used  to  calculate  the  amplitude 
fluctuations,  the  phase  fluctuations,  and  the  scattering  coefficient.  Theo- 
retical formulas  were  developed  for  the  fluctuations  and  the  scattering 
coefficient  tn  the  Fraunhofer  and  the  Fresnel  approximations  for  the  specular 
direction.^ 

The  theoretical  approach  of  Gulin  was  used  successfully  in  the 
FrasneL  approximation  to  predict  the  observed.-  amplitude  fluctuations  for 
value a of  the  Rayleigh  parameter  less  than  Q.6.  The  experimental  amplitude 
fluctuations  reached  a maximum  at  a Rayleigh  parameter  of  t.§  and  for  values 


between  if, § and  ’?.t  the  fluctuations  disclosed  a spread  of  % - §5  percent 
(relative  fc.e  the  average  amplitude). 

Gulin’s  theory  in  both  the  Fraunhofer  and  Fresnel  approximations 

agreed  with  the  experimental  phase  fluctuations  t@  within  6 percent  . 
(relative  t@  a radian)  for  values  of  the  Rayleigh  parameter  less  than  0.6. 

The  experimental  phase  fluctuations  varied  frets  12-  t©  £59  percent  (relative 
to  a radian)  over  a range  of  values  of  Rayleigh  parameters  from  9.2  to  7 «4. 

The  theoretical  development  in  the  Fraunhofer  approximation 
had  to  be  extended  to  the  Fresnel  approximation  In  order  t©  obtain  values  - 
for  the  scattering  coefficient  that  were  in  good  agreement  with  the  measured 


iv 


values  in  the  specular  direction.  As  the  value  of  the  Rayleigh  parameter 
increased  the  difference  between  the  theoretical  values  and  the  measured 
values  of  the  scattering  coefficient  increased  in  both  the  Fraunhofer  and 
Fresnel  approximations.  For  values  of  the  Rayleigh  parameter  near  1.8  the 
Fresnel  theory  predicted  values  within  2 d£  of  the  experimental  values  of 
the  scattering  coefficient.  When  th<*  Rayleigh  parameter  for  the  surface 
was  3*o  (7.2),  the  theoretical  value  for  the  scattering  coefficient  was 
3 dB  (10  dB)  smaller  than  the  experimental  value.  The  need  for  the  Fresnel 
theory  was  evident  in  tlxat  for  values  of  the  Rayleigh  parameter  near  7 .2 
the  Fraunhofer  theory  predicted  values  that  were  22  dB  smaller  than  the 
experimental  scattering  coefficient. 
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CHAPTER  I 
INTRODUCTION 


The  problem  of  the  scattering  of  waves  from  a rough  surface 
has  interested  physicists  ever  since  Rayleigh  (1896)  first  predicted 
the  scattered  field  from  a sinusoidal  surface.  His  solution  was 
limited  to  normal  incidence  and  was  applicable  only  when  the  wave- 
lengths of  the  surface  and  the  incident  radiation  were  much  larger 
than  the  surface  amplitude.  It  was  67  years  later  before  an  exact 
solution  for  the  problem  of  the  reflection  of  a plane  sound  wave  from 
a pressure-release  corrugated  surface  was  obtained  by  Uretsky  {1963). 

He  made  no  simplifying  assumptions  concerning  the  frequency  of  the 
incident  radiation  or  the  geometrical  parameters  of  the  sinusoidal 
surface.  Uretsky  (1965)  has  published  an  expanded  paper  in  which 
he  gives  more  details  regarding  his  solution  of  this  problem.  Barnard, 
Horton,  Miller,  and  Spitznogle  (1966)  have  given  experimental  results 
which  show  that  Uretsky  fs  theory  provides  a satisfactory  prediction 
of  scattered  sound  field  from  a pressure-release  sinusoidal  surface 
when  the  amplitude  of  the  sinusoid  is  comparable  to  the  wavelength 
of.  the  incident  radiation. 

The  number  of  papers  published'  in  the  area  of  underwater 
acoustics  in  the  last  15  years  has  shown  the  growing  interest  in 
understanding  the  scattering  of  underwater  sound  by  the  sea  bottom 
and  the  sea  surface.  The  topography  of  these  surfaces  is  neither 
periodic  nor  deterministic;  in  practice  they  can  only  be  described 
by  their  statistical  properties.  The  first  general  theory  which 
predicted  the  field  scattered  by  a perfectly  reflecting  statistically 
rough  surface,  the  irregularities  being  large  compared  with  the  wave- 
length, was  presented  by  Isakovich  (1952).  The  problem  was  solved 
in  the  Kirchhoff  approximation;  that  is,  it  was  assumed  that  the 
field  at  each  point  of  the  surface  could  be  represented  as  the  sum 
of  the  incident  wave  and  a wave  reflected  from  the  plane  tangent  to 
the  surface  at  the  given  point.  In  addition,  the  Fraunhofer  (farfield) 
approximation  was  used  throughout  the  paper.  The  rough  surface  that 


was  studied  by  Isakovich  was  isotropic  and  characterized  by  a Gaussian 
covariance  function. 

Eckart  (1953)  presented  a theory  which  predicted  the  field 
scattered  by  a pressure  release  or  rigid  surface  that- was  randomly 
rough  and  isotropic.  His  theory  for  a directional  source  leads  to 
useful  expressions  for  two  limiting  cases: 

(1)  when  the  wavelength  of  the  incident  radiation  is  much 

greater  than  the  largest  values  of  the  surface  relief,  and 

(2)  when  the  wavelength  is  much  less  than  the  surface  relief. 
The  Kirchhoff  and.  the  Fraunhofer  approximations  were  employed  through- 
out  the  paper.  Eckart  did  not  consider  any  specific  covariance 
functions. 

Apparently,  the  first  theory  in  the  Fresnel  approximation 
(as  opposed  to  the  Fraunhofer  approximation)  was  presented  by 
Feinstein  (1954).  The  method  of  physical  optics  (Kirchhoff  approxi- 
mation) was  used  to  predict  the  field  scattered  from  a randomly  rough 
surface  characterized  by  a Gaussian  covariance  function.  Also,  he 
employed  the  exponential  covariance  function  to  describe  the  field 
scattered  from  a one-dimensional  rough  surface . All  formulas  were 
restricted  to  the  specular  direction. 

Clay  (i960)  extended  the  theory  of  Eckart  (1953)  to  include 
an  omnidirectional  source.  A Gaussian  covariance  function  for  the 
rough  surface  was  assumed  for  the  calculations.  Considering  only  the 
long  wavelength  case,  Clay  compared  his  theoretical  rv,* suits  for  the 
amplitude  fluctuations  to  the  fluctuations  of  the  sound  reflected 
from  the  sea  surface  measured  by  Brown  and  Bicard  (i960).  The  numer- 
ical calculations  of  scattered  sound  had  the  same  dependence  on  tile 
source-receiver  separation  as  the  experimental  data. 

Formulas  for  the  amplitude  and  phase  fluctuations  of  a sound 
wave  reflected  from  a statistically  uneven  surface  were  derived  by 
Gulin  (IJoS).  His  theoretical  approach  was  similar  to  that  of  Sekart 
(1953)  except  that  all  calculations  were  made  in  the  Fresnel  approxi- 
mation and  were  confined  to  the  specular  direction.  Daly?  the  long 
wavelength  case  in  the  Kirchhoff  approximation  was  considered.  The 
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details  were  given  for  the  derivation  of  the  fluctuations  of  a 
spherical  wave  scattered  from  a rough  surface  characterized  by  a one- 
dimensioi  . quasi -harmonic  Gaussian  covariance  function,  The  details 
in  the  development  using  a two-dimensional  Gaussian  covariance  function 
were  omitted.-  The  theoretical,  results  for  the  amplitude  fluctuations 
obtained  with  the  one-dimensional  quasi-bermonie  Gaussian  covariance 
function  were  compared  with  the  results  of  an  experimental,  investi- 
gation made  by  Gulin  and  Malyshev  (1962)  of  th&;a^iitud^ 'fldotuatitip^^ 
of  sound  signals  reflected  from  the  undulating  surface  of  the  sea. 

The  experimental  amplitude  fluctuations  increase!  linearly  with  the  ■ • 
Rayleigh  parameter  up  to  values  of  the  latter  equal  to  approximately 
0.7  and  corresponded  to  the  theoretical  law  ever  this  range  of  Rayleigh 
parameters.  For  Rayleigh  parameters  between  0.7  and  3.5  there  was  a 
well  defined  alternation  of  experimental  amplitude  coefficient  of 
variation  between  maximum  and  minimum  values.  This  saturation  effect 
at  large  values  of  the  Rayleigh  parameter  could  not  possibly  be 
explaired  by  means  of  the  theoretical  analysis,  as  the  region  of 
applicability  of  the  analytical  formulas  was  limited  by  the  condition 
that  the  Rayleigh  parameter  must  be  less  than  one.  It  was  stated  by 
these  authors  that  the  conclusion  of  the  maxima  and  minima  in  the 
saturation  region  of  the  coefficient  of  variation  should  be  considered 
as  a tentative  result  and  in  need  of  further  verification.  No  experi- 
mental data  were  given  for  the  phase  fluctuation®. 

The  problem  of  the  scattering  of  waves  from  a rough ■ surface 
has  been  the  subject  of  at  least  one  monograph  written  by  Beckmann 
and  Sglaaiehifio  (1963).  This  hook  provides  m excellent  survey  of 
the  literature.  Beskaaan*®  approach  to  the  problem  ia  more  general 
than  that  of  Isakovieh  (1W&)  and  that  ef  gokort  (1953)  in  that  it  is 
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not  restricted  t©  tine  special  eases  ©f  a short  wavelength  or  a long 
wavelength  assumption.  Beckmans  employed  both  the  Kirehheff  and  the 
Fraunhofer  approximation  throughout  the  theoretical  analysis.  A 
rough  surface  characterised  by  .a  Gauss? an  ervari&aee  function  was 
investigated.  Later,  Beeh&aan  (lofi**)  discussed  the  scattering  of  a 
wave  feata  a tough  surface  characterised  by  an  exponential  emuiisee 


function.  Medwin  (1966)  modified  Beckmann's  theory  to  account  forthe-" 
error  inherent  in  the  common  assumption  of  uniform  insonification  over 
a fixed  area." 

Horton  and  Muir  ( 1967 ) have  employed  the  theory  of  Eckart 
(1953)  to  study  the  scattering  from  isotropically  rough  surfaces  whose 
covariance  'unctions  are  exponential,  Gaussian,  and  sinusoidal.  The  v- . . 
effect  of  the  radiation  pattern  of  the  source  upon  the  scattered 
intensity  was  discussed  by  inserting  an  analytical  expression  for  the 
transducer- pattern.  In  a companion  paper  Horton,  Mitchell,  and  Barnard 
(1967)  produced  surprising  experimental  evidence  of  the  validity  of  a 
modified  form  of  Eckart 's  theory  in  the  short  wavelength  case.  These 
experimental  r vudies  were  performed  on  a model  surface  whose  covariance 
.unction  was  exponential . The  theory  was  developed  in  the  diirchhoff 
and  the  Fraunhofer  approximations. 

Of  the  references  quoted  here  only  two,  Feinstein  (195*0  and 
Gulin  (ig6i0,  have  studied  the  problem  of  scattering  from  a rough 

N.' 

surface  in  the  I esnel  approximation.  Neither  have  compared  their 
results  with. -the  correspor ling  solutions  in  the  Fraunhofer,  approxi- 
mation. One  migfct  question  whether  or  not  the  extension  of  the 
mathematical  analysis  of  the  problem  from  the  Fraunhofer  to  the 
Fresnel  approximation  world  give  better  agreement  between  the  theoreti- 
cal jmd  -experimental  re- ‘fits.  No  one  to  tne  author’s  knowledge  has. 
presented  a solution  in  the  Fresnel  approximation  te  the  problem  of 
scattering'  from  a randomly  rough  surface  who m statistical  properties 
are  f'iven  by  an  exponential  covariance  function . One  is  further 
motivated  to  investigate  the  scattering  of  acoustic  waves  from  this 
type  of  surface  since  there  are  available  model  surfaces  that  are 
characterised  by  an  exponential  covariance  function.  Those  model  sur- 
faces whose  statistical  properties  were  analysed  by  Horton,  Hoffman, 
and  Heakins  (1962}  give  tu  excellent  mesas  of  experimental  verification 
of  the  theoretical  results. 

The  plan  of  study  will  be  the  following.  In  Chapter  XI  - 
general  expressions  are  derived  for  the  amplitude  and  phase  fluctuations 
of  an  acoustic  wave  scattered  in  the  specular  direction  i Yob  a randomly 


rough  surface.  The  theoretical  approach  is  that  of  Gulin  (1962).  In 
the  same  chapter  general  expressions  are  obtained  for  the  scattering 
coefficient  derived  in  a manner  similar  to  that  presented  by  Beckmann 
{1963}.  Throughout  Chapter  II  there  is  a parallel  presentation  of  the 
integral  expressions  obtained  in  the  Fraunhofer  and  Fresnel  approxi- 
mation. The  expressions  for  the  amplitude  and  phase  fluctuations  will 
be  valid  only  for  Rayleigh  parameters  less  than  one  whereas  the 
^formulas  for  the  scattering  coefficient  should  be  correct  for  any 
value  of  the. Rayleigh  parameter.  The  calculations  in  this  initial 


r?ftudy  shall  be  confined  to  ^1i%^sp'eculai^dipe<|tion  where  the  mathe- 


iiaatical  analysis  is  less  cumbersome. 

\ , In  Chapter  III  the  problem  is  specialized  to  the  case  where 

•» 

the  statistical  properties  of  the  rough  surface  are  described  by  an 
exponential  covariance  function.  The  integral  expressions  occurring 
in  Chapter  II  are  solved  giving  formulas  for  the  amplitude  fluctuations 
phase  fluctuations,  and  scattering  coefficient  that  can  be  used  to 
describe  the  experimental  results. 

The  statistical  properties  of  the  model  surfaces  used  in  the 
experimental,  work  are  discussed  in  Chapter  IV.  This  section  presents 
a list  of  the  experimental  equipment  and  gives  an  explanation  of  the 
experimental  procedures  used  to  measure  the  fluctuations  and  the 
scattered  intensity. 

The  experimental  data  are  compared  with  the  theoretical 
formulas  in  Chapter  V,  a general  discussion  is  presented  in  Chapter  VI, 
and  the  conclusions  of  this  s*udy  are  stated  in  Chapter  VII. 

The  theory  presented  in  Chapters  II  and  HI  and  in  Appendixes 


A and  C can  best  be  followed  by  referring  to  the  block  diagrams  in 
Figs,  1 and  g. 
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CHAPTER  II 

THEORETICAL  APPROACH 

General  expressions  are  derived  for  the  amplitude  and  phase 
fluctuations,  and  for  the  scattering  coefficient  in  the  specular  direction 
for  an  acoustic  wave  scattered  from  a randomly  rough  surface.  The  theo-  . 
retical  approach  that  is  used  is  the  method  presented  by  E.  P.  Gulin,  (1962} . 
Gulin's  notation  will  be  maintained  in  this  paper. 

A.  The  Pressure  Field 

The  coordinates  of  the  source  (Q,0,z^)  and  the  receiver  (L,Q,Zg) 
are  illustrated  in  Fig.  3.  ^ is  the  grazing  angle.  ^ is  the  receiving 

angle  measured  with  respect  to  the  x-axis.  The  axis  of  the  directional 
source  is  directed  toward  the  point  on  the  surface  labeled  O' . The 
x-y  plane  is  oriented  so  that  the  average  height  of  the  interface  above 
this  plane  is  zero.  is  tue  distance  from  the  source  to  the  point  O' 
which  is  located  on  the  x-y  plane.  is  the  distance  from  the  source  to 
any  arbitrary  point  on  the  randomly  rough  surface.  is  the  distance 

from  the  source  to  a point  on  the  x-y  plane  directly  above  or  below  the 
arbitrary  point  on  the  scattering  surface.  The  distances  R^»  arMi  Kg 
to  the  receiver  are  defined  in  the  same  way. 

The  total  pressure  field  p of  a sound  wave  reflected  from  an 
uneven  surface  can  be  written  in  terns  of  the  Helmholtz  Integral,  Baker 
and  Copsoa  (1939)» 


p * h jj  h k 


3a  45 


where  n Is  the  outward  normal  to  the  surface  (in  a direction  away  from 
source  and  receiver).  The  wave'  incident  on  the  surface  is  assumed  to  be 
a spericai  wave  given  by 
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FIGURE  3 

COORDINATE  SYSTEM 


The  field  at  the  surface  S is  specified  after  the  example  of 
Eckart  (1953)  for  the  special  case  of  a "pressure  release"  surface.  The 
boundary  conditions  at  the  surface  S are 


o'  + p = 0 , and 

o 


(3) 

W 


The  application  of  these  boundary  conditions  corresponds  to  the  Kirchhoff 
approximation  in  the  problem  of  sound  scattering  by  rough  surfaces.  The 
applicability  of  +his  approximation  has  been  formulated  by  Brekhovskikh 
(1952)  and  Isakovich  (1952),  and  has  been  discussed  by  Beckmann  (l^^v'.-,;3 
Substituting  Eqs.  (2),  (3),  and  (4)  into  Eq.  (l)  one  obtains 
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p c 


Id R.  dRj 
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ik(Rx+Rg) 
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where  >.  is  fete  wavelength  of  the  incident  radiation.  If  it  is  assumed 


Iv^  » X , and 

it,  » \ 


(6) 


the  second  group  of  toms  can  be  neglected . which  reduces  Eq.  0)  to 


ik 


* \** 

4 


du  1 


dS 


(T) 


in  the  case  of  1 rregtiiarit les  ©?  sufficiently  small  slope' 
differentiation  along  the  normal  can  be  approximately  replaced  by 
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differentiation  along  the  negative  z-axis,  and  integration  over  the 
rough  surface  by  integration  over  the  x-y  plane/  Thus,  Eq.  (7) 
becomes 


+ 

3T  Fi 


ik(VR  ) 
e x * dxdy 


(8) 


Introducing  the  new  terns  z£,  F which  are  depicted  in 

Fig.  4,  one  can  show  that 

asi_  “i  < * 

■3s  Rx  ' 3s“"“Rg  ■* 

Eq.  (8)  becomes 


ik(R 

e dxdy 


(9) 


F to  the  deviation  of  the  surface  relief  from  the  average  height.  If  it 
is  as siued  that  |P|  is  much  less  than  the  minimum  of  and  s^,-  then 
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%ceteaa«  (iPJ)  inserts'  tM  “local'1  seatterittg  of  the  ncraal 

derivative  at  tins  1®  his  theoretical  approach,  this  method 

introduces  the  tests 
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which  reduces  to  on*  la  the  specular  direction  = % ), 


FIGURE  4 

RANDOM  QUANTITIES  z's'P 


Defining  (see  Fig.  5) 
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one  can  write  as 


Similarly 


*2  “ ■ X^2  + yC  * <v:t)3 


and  Rg.  can  tie  written  as 
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Since  P « z1  and  P « Zg,  the  second  tern  in  the  square  brackets  can  be 
neglected,  ‘fhe  pressure  becomes 


ik  rr.^1 

p”^  jj  r ‘ k si 
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- Z A “ikP  IP*  + 
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If  the  acoustic  wavelength,  is  greater  th>»n  the  surface  relief 
F (x,y),  or  if  the  grasiug  angle  is  small,  so  that  the  following  condition  ;; 
holds  ' 


ri  zs 
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21 
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This- approximation  separates  Sq.  (l4)  into  a aonfluetuating  eoaspoaent 


"91  1 


dsdy.  , . 


aM  a scattered,  fluctuating  eamjp&s&t 
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She  f®egsdf'®  p reflected  frees  a plane  "p-essurc  release*  surface  is 
evaluated  for  the.  specular  directies  in  ApeeMi*  -C. 


B.  Amplitude:  and  Phase  Fluctuations  in  the  Fraunhofer  Approximation 

The  transition  from  the  expressions  for  the  sound  field  to 
expressions  for  the  amplitude  and  phe.se  fluctuations  can  be  carried  out 
by  using  a method  utilized  by  Chernov  (i960).  The  total  pressure  given 
by  Eqs.  (IT)  and  (l8) 


can  be  written  as 


P=  P0  + P-l 


i®, 

Ae1$  = AQe  0 + A^e.  1 


i(®-$  ) A if®..-®  ) 
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= l + f 6 
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Assuming  that 


one  can  set 


A1«Aq  , 


A = A + SA  , and 
o 7 


® = ® + 6® 
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These  approximations,  with  the  use  of  the  logarithm  and  in  (l  + e)  a € 
when  je|  « 1,  reduce-  Eq,.  (19)  to 


A,  if®,-®  ) 
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Since  the  right  side  of  this  equality  is  Eq.  (l8)  divided  by  Aa  , one 


has  for  the  variation  in  the  amplitude 
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k(Rj.  + Rg)  - ®Q  F(x,y)dxdy  , (20) 


and  for  the  variation  in  phase 


r 

l2i  + “21 

jJ 

RiR2 

K -SaJ 

jj  m re  w 

In  the  Fraunhofer  approximation  (Appendix  A,  Eq.  (A- 20)  for  t 


specular  direction 


Ri  + R*  * Ri0  + Hgo 
so  that  the  variation  in  amplitude  is 


k2  sin2  t cos[k(R-^Q  + ^20^  ~ ^o] 

*Ao*l(p20 


J J F(x]_,y).^1dy  . (22), 


where  (Appendix  A,  Eq.  (A-l) 


X1  = X ’ ^ + z2 


The  variation  in  phase  is 
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k sin 
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nAoRl0R20 


J J FU-^y^dy 


, w (r  and  (C-M  in  Appendix  C.  These  values 
A and  $ are  given  by  Eqs.  (C-3)  and  . I Vir' 

reduce  Eqs.  (22)  and  (23)  to 
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where  A^,  A^,  and  l are  the  dimensions  of  the  active  scattering  area  defined 
in  Appendix  B.  5$  is  zero  only  if  the  integration  extends  over  a signifi- 
cant area  of  the  surface.  If  the  active  scattering  region  is  small  compared 
to  the  correlation  distance,  8$  is  likely  to  differ  from  zero  for  each 
single  measurement. 

Since  F(x^,y)  is  a random  variable,  the  quantities  of  interest 
are  the  ensemble  averages.  It  is  assumed  that  the  surface  relief  is 
distributed  so  that 


(26) 


Since  the  variation  in  amplitude  is  identically  zero,  one  can  write 


(27) 


without  reference  to  the  distribution  of  the  surface  relief. 

The  variance  of  a random  variable  v is  defined  as  the  average 
square  of  the  deviation  from  the  mean  value 


D(v)  = ^(v  - <v»2)  = - (v)2  . (28) 

The  square  root  of  the  variance, 

a = /D(vl  . , (29) 

is  called  the  standard  deviation.  Thus,  utilising  Eq,  (27),  one  can  write 
for  the  relative  standard  deviation  of  the  amplitude  the  expression 


which  hereafter  will  be  called  the  amplitude  fluctuations.  In  the  same  way 


will  represent  the  phase  fluctuations. 


-It 


Eq.  (24)  indicates  that  the  amplitude  fluctuations  in  the 
specular  direction  are  zero  in  the  Fraunhofer  field  zone.  To  determine 
the  phase  fluctuations  in  the  specular  direction  in  the  Fraunhofer 
approximation  Eq.  (25)  is  squared  and  averaged  to  give 
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(52) 


where  ^(x^y^x^y'))  is  defined  as  the  spatial  covariance  function 
of  the  surface  displacements.  A complete  solution  of  the  problem  will  be 
delayed  until  Chapter  III  where  by  Eq.  (56)  an  exponential  convariance 
function  is  specified. 

C.  Amplitude  ana  Phase- Fluctuations  in  the  Fresnel  Approximation 

The  Frssnel  approximation  for  the  expansion  of  in  the 

specular  direction  is  given  by  Kq.  (A- 19)  in  Appendix  A, 


where 
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Substitution  of  Eq.  (53)  into  Eq.  (20)  gives  for  the  variation  in  amplitude 
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Aq  and  $q  are  given  by  Eqs.  (C-10)  and  (C-ll)  in  Appendix  C.  Eqs.  (33) 
and  (21)  give  for  the  variation  in  phase 
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F(x1,.y)dx1dy 


(37) 


Assuming  that  the  surface  relief  is  distributed  such  that 
the  ensemble  average  of  Eqs.  (35)  and  (37)  gives 

= 0 , and 


(39) 


It  is  of  interest  to  noint  out  that  in  the  Fraunhofer  approximation  SA/A_ 

s- 

is  identically  aero,  whereas  in  the  Fresnel  case  it  is  only  the  average 
value  («/y  that  vanishes.  '• 

Squaring  and  averaging  Eq.  (35)  and  utilising  Eq.  (38)  one 
obtains  for  the  square  of  the  amplitude  fluctuations 
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Eqs.  (37)  and  (59)  give  for  the  square  of  the  phase  fluctuations 
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The  development  of  Eqs.  (40)  and  (4l)  will  be  continued,  in  Eqs.  (67)  and 
(75)  of  Chapter  III,  respectively. 

D.  Scattering  Coefficient  in  the  Fraunhofer  Approximation 

The  Fraunhofer  approximation  in  the  specular  direction  is  the 
substitution  from  Eq.  (A-2Q)  in  Appendix  A of 

*1  * ^ " ®10  + ll2Q 

into  the  exponential  of  Eq.  (14).  With  ^ sin  ^ ands^  = oin  it 

Eq.  (14)  beeottea 


- ik  (sin  #)e  ^ f f .-21 
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The  intensity  e-f  the  scattered  wave  is  defined  as 

k 


where 


a la  the  density, 

v is  the  velocity  of  the  aedius  through  which  the  wave  is 
p is  the  cecspiex.  conjugate  of  the-  pressure  p. 


•'x  ■. 


e: 


The  product  pv  is  called  the  specific  acoustic  impedance  of  the  medium. 
Forming  the  product  pp  from  Eq.  (42)  and  taking  the  ensemble  average,  oueb 
obtains 


(pp*)  = 
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It  is  assumed  that  the  surface  relief  of  the  rough  surface  is 
characterized  by  a bivariate  Gaussian  distribution  function,  Eckart  (1955) 
and  Horton  and  Muir  ( 19^7 ) , defined  by 


where 


C * ^F(x1,y)F(x^y* 
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is  the  spatial  eovarianee  funetion  of  the  surface  displacements,  and 


he  shown  that  with  this  function  one  obtains 


s e 
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Eq.  (46)  reduces  Eq.  (43)  to 

2 2 -4k2  /p2\sin2y 

k sin  \ o / 

1 - 2_  2 „ 2 e 
4n:  HgQ 

(47) 
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The  scattering  coefficient  ag  is  defined  as  the  ratio  of  the 
intensity  scattered  from  a rough  surface  to  the  intensity  reflected  from 
a plane  surface. 
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The  tern  p p , which  ia  proportional  to  the  intensity  reflected  fttsi  a 

Q 0 

plane  surface,  is  given  in  Appendix  C by  Eq.  (C-S)» 
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where  4|,  and  I are  the  dimensions  of  the  active  scattering  region 
diseased  in  Appendix  3.  Thus,  the  scattering  cue  f fie  lent  in  the 
Fraunhofer  approximation  in  the  specular  direction  is  given  fey 
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The  intensity  reflected  from  a.  plane  surface  in  the  Fresnel  approximation 
is  proportional  to  the  expression  given  by  Eq,.  (C-12)  in  Appendix  C, 
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The  Fresnel  Integrals , § and  C,  art*  defined  by  Kq.  (C-9>  in  Appendix  0. 

fhe  scattering  coefficient  is  the  ratio  of  Sis.  (52)  and  (5S)» 
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Kq.  (55)  will  be  carried  forward  in  Eq.  (90)  of  Chapter  III 
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CHAPTER  III 

THEORY  USING  AN  EXPONENTIAL  COVARIANCE  FUNCTION 


The  statistical  properties  of  the  model  surfaces,  which  shall  he 
discussed  in  Chapter  IV,  suggest  that  an  exponential  function  best  describes 
the  spatial  covariance  of  the  rough  surfaces  employed  in  the  experimental 
study.  The  exponential  spatial  covariance  function  of  the  surface  dis- 
placements is  given  by 


c=  = (?iy 


-l|<Vxi 


)2+{y-y*)^ 


where  a is  the  spatial  correlation  interval  of  the  surface  displacements 

and  is  the  mean-square  amplitude  of  the  surface  irregularities. ■ The 

distance  a,  for  which  C will  drop  to  the  value  e”1,  has  the  physical  sig- 
nificance that  it  is  a measurement  of  the  relative  distance  between  the  hills 
and  valleys  on  the  randomly  rough  surface.  The  square  ro't  of  (f2^  has  the- 
physieai  meaning  that  *t  represents  the  average  heights,  of  the  hill3  and  the 
average  depth  of  the  valleys. 

A.  AmpUtudf*  nnl  Flvu-e  Fluctuations  in  the  Fraunhofer  Approximation 

The  amplitude  fluctuations  in  the  Fraunhofer  approximation  are  aero 
irrespective  of  the  form  ©f  the  spatial  covayiar.se  function  as  shewn  by 
Eq.  (ik)  in  Shifter  II. 

The  phase  fluctuations  in  the  Fraunheftr  approximation  are 
determined  by  the  substitution  ©f  Sq.  (56)  into  Eq.  (5§)  te  give 
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CHAPTER  III 

THEORY  USING  AN  EXPONENTIAL  COVARIANCE  FUNCTION 

The  statistical  properties  of  the  model  surfaces,  which  shall  he 
discussed  in  Chapter  IV,  suggest  that  an  exponential  function  best  describes 
the  spatial  covariance  of  the  rough  surfaces  employed  in  the  experimental 
study.  The  exponential  spatial  covariance  function  of  the  surface  dis- 
placements is  given  by 


/ \ / 2\  -ak-*i>2*<y^>2]  ^ 

C = <F(xI>y)F(x[,y,y  = Yo/e  *-  , (56) 

where  a is  the  spatial  correlation  interval  of  the  surface  displacements 

and  Yq)  the  mean-square  amplitude  of  the  surface  irregularities.  The 
distance  a,  for  which  C will  drop  to  the  value  e~\  has  the  physical  sig- 
nificance that  it  is  a measurement  of  the  relative  distance  between  the  hills 
and  valleys  on  the  randomly  rough  surface.  The  square  r ,t  of  (?~^  has  the 
physical  meaning  that  {t  represents  the  average  heights  of  the  bills  and  the 
average  depth  of  the  valleys. 

A.  Amplitude  an .1  Phase  Fluctuations  in  the  Fraunhofer  Approximation 

The  amplitude  fluctuations  in  the  Fraunhofer  approximation  are  zero 
irrespective  the  form  of  the  spatial  covariance  function  as  shewn  by 
Kq.  (24)  in  Chapter  II. 

The  phase  fluctuations  in  the  Fraunhofer  approximation  are 
determined  by  the  substitution  of  Eq.  ) into  Eq.  (Jt)  te  give 
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Since  and  Ag  are  not  equal,  it  is  convenient  to  make  the  preliminary 
change  of  variables  (see  Fig.  5) 

A2  “ A1 

..  , . u=  ^ g—  > 


. A2  “ *1 
u = *1 


and  introduce 


W-. 


With  these  changes  Eq.  (57)  becomes 
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fhe  approximate  sise  of  thf  active  scattering  reiien  is  by  Si  oa  defined 
in  Appendix  B. 

If  one  Introduces  relative  geordin&teg  as  shews  in  Fig.  5 


defined  by 
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■ asi  eente'r  of  gravity  coordinates  defined  by 
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Eq:  '60  )vbecomes ' 
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(63) 


: The.  integrals' in  Eq.  (63)  can  not  be  computed  in  open  form. 

Further  discussion  shall  be  restricted  to  the  case  when  the  dimensions  of 
the  region  essential  to  integration  over  | and  n are  considerably  less  than 
those  of  the  irradiated  part  of  the  surface.:  That  is,  it  is  assumed  that 

a « 2A,  2£  (6h) 

In  keeping  with  this  approximation  one  may  extend  the  range  of  integration 
over  i and  tj  to  .infinity,  so  that  Eq.  (63)  becomes 
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The  integration  may  be  performed  readily  by  changing  to  polar  coordinates, 
and  one  finds  for  Eq.  (65)  ' 


a.  - 
0 


(66) 


this  equation  shows  that  the  phase  fluctuations  in  the  Fraunhofer 
approximation  are  proportional  to  the  -Rayleigh  parameter,  2k  j (V 2 ^ sin  ^ 

(for  Rayleigh  parameters  much  less  than  one),  and  to  the  spatial  correlation 
interval,  &»  The  fluctuations  are  inversely  p r op 0 r t i on a 1 to  the  square  root 
of  the  approximate  insonified  area,  21  ( + A^) 


Substituting  Eq.  (56)  into  Eq.  (**0)  of  Chapter  II,  one  finds  for 
the  square  of  the  amplitude  fluctuations 
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The  details  in  the  integration  of  J are  given  in  Appendix  D.  The  result  is 
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where  C and  S are  the  Fresnel  integrals  defined  by  Eq.  (C-9)  in  Appendix  C. 
V^y.-  A discussion  of  the  term  is  given  in  Appendix  E.  The 

v investigation -.  shows  that 


V J « J 
y2  1 


(71) 


so  that  the  amplitude  fluctuations  as  given  by  Eq.  (6?)  can  be  approximated 
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Substituting- Eq.  (56)  into  Eq.  ;(4l)  of.  Chapter  II,.  one  obtains 
. for  the  phase  fluctuations 


/{(M)2)  = /j ~Jr 


(73) 


where  and  are  given  by  Eqs.  (68)  and  (69).  Since  J g « J^,  the 
amplitude  and  phase  fluctuations  are  approximately  equal  for  small  Eayleigh 
parameters  so  that 
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Thus,  the  Fraunhofer  approximation  shows  the  phase  fluctuations 
to  be  significantly  larger  than  the  amplitude  fluctuations  (which  are  zero) 
whereas  the  Fresnel  approximation  shows  the  amplitude  and  phase  fluctuations 
to  be  approximatley  of  the  same  magnitude.  - 

C.  Scattering  Coefficient  in  the  Fraunhofer  Approximation 

Substituting  Eq.  (56)  into.  Eq,  (^9)  of  Chapter  II,  one  obtains 
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g a 4k2  sin“  i|r  (76) 

and  refer  to  g as  the  Beckmann  parameter  (which  in  the  specular  direction  is 
the  square  of  the  Rayleigh  parameter). 

The  integrand  can  be  written  as 
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Eq,  (75)  becomes 
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Since  the  integrals  occurring  in  Eq.  (78)  are  the  same  as  those 
occurring  in  Eq.  (57)  for  the  phase  fluctuations,  one  can  write  immediately 
from  the  square  of  Eq.  (66), 
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for  the  scattering  coefficient  in  the  specular  direction  in  the  Fraunhofer 
approximation.  '„■  ■ - - 


Scattering  Coefficient  in  the  Fresnel  Annroximation 


Substituting  Eq.  (56)  into  Eq.  (55)  of  Chapter  II,  one  obtains 
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Using  the  infinite  summation  defined  by  Eq.  f 77) » one  obtains 
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Integrating  the  first  set  of  integrals  as  shown  in  Appendix  C 
and  separating  the  second  set  of  integrals  into  their  real  and  imaginary 
parts,  one  obtains 
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The  integration  of  the  second  group  of  integrals  containing  the  sine 
function  can  be  accomplished  by  following  the  procedure  given  for  the 
integration  of  the  J1  tens  in  Appendix  D.  Tfee  parallel  development 
produces  the  integral  (see  £4.  (1D~3)  in  Appendix  D) 
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which  can  be  reduced  to 
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From  Gradshteyn  and  Byahik  (integral  3.715-15,  page  402),  one  obtains 

1=0.  Thus,  the  expression  for  the  coefficient  o becomes 
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The  integrals  appearing  in  Kq.  (82)  are  the  same  as  the  Integrals  appearing 
in  (0).  thus,  the  scattering  coefficient  calculated  in  the  Fresnel 
apprexisati&n  can  be  written  as 


in 


ft  M g 

®s 


l.f  ) j£j(S) 

g / ta.  1 o 
hpI 


where  jL  (§)  is  given  by  E%.  (70). 
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CHAPTER  IV 

EXPERIMENTAL  APPARATUS  AMD  PROCEDURES 

It  is  seen  from  the  theory  developed  in  Chapter  II  that  one  must 
know  the  statistical  properties  of  the  rough  surface  in  order  that  the 
predictions  of  the  theory  can  be  compared  with  measurements  on  the  rough 
surface.  In  particular,  if  F(x^,y)  is  the  elevation  of  the  reflecting 

surface,  one  must  know  the  autocovariance  <F(x1y)F(*j/y!  as  shown  by 

Eqs.  (52),  (40),  (4l),  (49),  and  (55),  and  the  bivariate  distribution 
function  as  illustrated  by  Eq.  (44). 

The  problem  of  selecting  a topographic  surface  to  be  modeled 
has  been  discussed  by  Horton,  Mitchell,  and  Barnard  (1967).  The  model 
surface  was  generated  from  maps  of  aeromagnetic  intensity  over  a portion 
of  the  Canadian  Shield.  Before  these  maps  were  used  to  construct  the 
model,  the  two-dimensional  autocovariance  functions  of  these  maps  were 
studied  by  Horton,  Hempkins,  and  Hoffman  (1964).  Details  of  the 
mathematical  analysis  used  in  the  investigation  of  such  maps  are  given 
by  Horton,  Hofftaan,  and  Hempkins  (1962). 

A summary  of  tiie  results  obtained  in  these  investigations  will 

a 

be  given.  The  Canadian  maps  were  contoured  in  gammas  (ly  >*  10  Qe). 

The  vertical  relief  of  the  model  surface  was  obtained  from  this 
aeromagnetic  dimension  using  the  scale  lQy  * 1/52  in.  The  factor 
1 mile  a 1 in.  was  used  to  convert  the  area  of  the  aeromagnetic  maps 
into  model  surfaces  of  site  52  in.  by  32  inches . The  theory  in 
Chapters  II  and  III  Is  based  on  the  assumption  that  the  reflecting  surface 
is  pressure  release  (see  Eq,  (3),  and  the  models  were  constructed  from  a 
material  of  small  acoustic  impedance  compared  with  the  acoustic  impedance 
of  tiie  water  in  which  the  surfaces  are  submerged.  The  model  was  made  of 
low-density  expanded  polystyrene  by  the  Construction  Service  Company  of 
Fort  Worth,  Texas.  A square  grid  of  sample  points  drawn  1/2  in.  apart  on 
the  central  one- fourth  of  the  map  surface  was  used  to  compute  the  auto- 
covariance  function.  This  square  grid  gave  106$  values  of  the  elevation. 


The  autocovariance  function  was  computed  from  a formula  given  by  Horton, 
Hoffman,  and  Hempkins  (1962) 

N-s  M-t 

C(M>  - ’01— I I F<1-»F<1  + ». » w 

i=l  J= 1 

where 

K = M = 33  is  the  number  of  samples  alone  each  l6  in.  grid  line, 

F is  the  surface  relief, 

s and  t are  units  of  shift  along  the  grid,  and 
C is  the  autocovariance  function. 

Profiles  of  the  auto covariance  function  along  the  e and  t directions  are 

plotted  in  Fig.  6.  The  profiles  show  that  the  surface  is  possibly 

isotropic  for  values  of  s and  t less  than  one  inch.  The  mean-square 
/ o\  ~ " p 

surface  relief,  (F0  ) » is  0.132  in.  A curve  that  will  pass  through 
the  three  points  nearest  the-  origin  and  approximately  average  the  other 
anisotropic  points  for  values  of  s and  t greater  than  one  inch  is  the 
exponential  autoeovarianee  function  given  by 

r 

C(r)  - 0.132  , (85) 

where  2.55  in.  1©  the  correlation  length  defined  as  the  distance  over 
which  the  covariance  function  is  reduced  by  l/e,  Eq,  (85)  corresponds 
to  Eq.  (36)  in  Chapter  III, 
where 

* 0.15?  ift.& 


a « t,55  in. 


Although  th  properties  of  the  model  surface  do  not  satisfy  the  assumption 
©f  Isotropy  on  which  the  theory  in  Chapters'  II  and  III  is  based,  the 
agreement  is  close  enough  to  justify  a caparison  between  theory  and 
experiment.  Since  no  one  has  measured  experimentally  the  bivariate 


_39. 


probability  density  for  a sample  of  topography,  it  is  not  possible  to 
justifv  directly  the  choice  given  by  Eq.  (44)  at  the  present  time. 
However,  Horton,  Mitchell,  and  Barnard  (1967)  give  convincing  reasons 
why  the  bivariate  Gaussian  distribution  may  be  used.  The  reader  is 
referred  to  a thesis  by  Mitchell  {1966)  for  further  details  concerning 
the  model  surface. 

Three  other  model  surfaces  described  statistically  by  the  same 
exponential  covariance  function  were  constructed.  The  statistical 
properties  of  the  four  model  surfaces  are  given  in  Table  X.  The 


ratio 


was  calculated  for  a frequency  of  95-800  kHz 


(wavelength  \ = 0.608  in.). 

The  experiments  were  performed  with  the  model  surface  submerged 
in  a water- filled  tank  with  dimensions  15  ft  wide,  60  ft  long,  and 
12  ft  deep.  The  surface  was  suspended  vertically  with  its  center 
approximately  46  in.  below  the  waterline.  An  x,y,z  mechanical  positioning 
system  mounted  on  the  top  of  the  tank  was  used  to  position  the  source 
and  receiver  in  the  experiment.  Figure  7 shows  the  laboratory  tank, 
positioning  system,  and  electronic  equipment  used  in  the  study.  The 
source  and  receiver  were  mounted  on  two  identical  columns  with  the  move- 
ment in  the  x,y,s  directions  read  from  a counter  (n  milli-inehes  per 
smallest  division  or  518.310  divisions  per  inch  on  the  counter).  Also, 
the  mechanism  allows  the  projector  and  receiver  to  be  rotated  about  the 
vertical  axis  in  increments  of  0.1  deg. 

The  source  was  a transducer  ttiat  has  an  active  face  diameter 
of  12.3  em  and  a resonance  frequency  of  94.0  kite.  The  directivity  pattern 
of  the  piston  transducer  is  shown  in  Fig.  6.  The  Ueanwidth  is  approxi- 
mately § deg  measured  between  points  where  the  radiated  intensity  falls 
to  one-half  the  intensity  along  the  forward  axis  (-5  dfi  points). 

The  receiver  was  a small  cylinder  of  outside  diameter  0,082  in., 
inside  diameter  0.05?  in.,  and  length  0,123  inch.  Its  directivity 
pattern  is  approximately  omnidirectional  as  shown  by  Fig.  9*  The  same 


side  of  the  probe  was  directed  along  the  slant-line  distance  to  the  model 
surfaces  throughout  the  experiment. 


*. 


TABLE  I 

STATISTICAL  PROPERTIES  OF  THE  MODEL  SURFACES 
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TABLE  I 

STATISTICAL  PROPERTIES  OF  THE  MODEL  SURFACES 


FIGURE 


S vi 


SOURCE  AND  RECEIVER  SEPARATED  BY  10  ft 
INPUT  VOLTAGE  IS  V 
PULSEtelOTM;  S40.  mwc 
RESONANCE  FREQUENCY:  94.0  kNs 


FIGURE  8 

DIRECTIVITY  PATTERN 
TRANJOUCER  DRL  No.  II602S?) 
FREQUENCY:  95.800  hHs 


180° 


SOURCE  AND  RECEIVER  SEPARATED  BY  10  ft 
INPUT  VOLTAGE:  IS  V 
PULSEWIDTH:  640  /xsec 
RESONANCE  FREQUENCY:  900  kHz 


FIGURE  9 

DIRECTIVITY  PATTERN  OF  RECEIVER 

FREQUENCY*  95.800  kHz 


Figure  10  is  a block  diagram  of  the  arrangement  of  the 
electronic  and  acoustic  equipment.  The  continuous  oscillator  signal  is 
gated  into  signals  of  640  psec  duration  at  a rate  of  28  times  per  second. 
These  pulses  are  95  cm  * 37  in.  long  in  the  water  and  contain  61  cycles 
at  95.800  kHz.  The  pulsed  signals -are  used  as  the  input  of  the  power 
amplifier  that  drives  the  projecting  transducer.  The  received  signal 
scattered  from  the  pressure  release  model  surface  is  amplified  and  is 
displayed  on  the  oscilloscope.  The  peak- to- peak  voltage  of  the  received 
signal  is  recorded.  At  the  same  instant  the  relative  phase  of  the 
received  signal  is  measured  by  recording  the  time  delay  of  the  scattered 
pulse  relative  to  a continuous  wave  from  the  oscillator,  The  counter 
allows  the  experimenter  to  maintain  the  frequency  to  within  oue  o'j  ole 
throughout  the  data  collecting. 

All  measurements  were  made  with  the  projector  and  receiver 
in  the  specular  direction.  With  the  projector  placed  60. 0 in.  from  the 
origin  on  the  rough  surface  and  with  the  receiver  placed  60,0  in.  frcan 
the  same  origin  (the  origin  lies  on  an  imaginary  plane  that  passes 
through  the  average  height  of  the  rough  surface),  the  grazing  angle  was 
varied  from  80  deg  to  6 deg  in  increments  of  2 deg.  Figure  11  shows 
samples  of  the  oscilloscope  pictures  for  grazing  angles  of  50  and  80  deg. 

The  center  of  the  received  pulse  is  recorded  since  this 
represents  a scattered  pressure  due  to  the  total  insonified  area.  The 
pulse  to  the  left  of  the  largest  signal  in  Fig.  11a  for  a grazing  angle 
of  80  deg  is  attributed  to  the  side  lobes  of  the  projector. 

If  the  source  and  the  receiver  are  displaced  as  a rigid  body 
parallel  to  the  scattering  surface,  oue  would  obtain  different  values 
for  the  amplitude  and  relative  phase  as  the  dominant  scattered  energy 
returns  from  different  parts  of  the  rough  surface.  Consequently,  the 
sequence  of  values  for  the  amplitude  and  relative  phase  constitute 
samples  of  a stochastic  variable.  Therefore.,  one  is  concerned  not  with 
specific  values  but  rather  with  average  values.  In  the  process. of 
collecting  the  data,  ten  different  values  of  the  amplitude  and  relative 
phase  were  obtained  by  insonifying  tea  different  regions  of  the  model 
surface.  To  inaoniiV  the  different  regions,  the  projector  and  receiver 
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EXPERIMENTAL  EQUIPMENT 


were  moved  as  a rigid  body  in  the  vertical  direction  to  ten  different 
positions  in  increments  of  one  inch.  Figure  12  shows  three  of  the  ten 
different  positions  of  the  area  that  the  projector  insonifies  (numbering  1, 
6,  and  10) » The  ellipses  indicate  the  approximate  size  of  the  insonified 
area  (-3  dB)  for  a grazing  angle  of  40  deg.  The  active  scattering  area 
for  each  of  the  ten  data  runs  is  confined  to  a rectangular  region  14.8  in. 
by  18.4  inches. 

The  amplitude  fluctuations  were  measured  for  any  one  grazing 
angle  with  the  receiver  confined  to  the  specular  direction  by  recording 
the  peak-to-peak  voltage  of  the  received  signal  for  the  ten  different 
positions  of  the  transmitter  and  receiver.  The  average  amplitude 
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was  used  to  calculate  the  variance 
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from  which  the  standard  deviation 


was  obtained,  The  relative  standard  deviation  (the  amplitude  fluctuations) 
io  the  ratio  of  the  standard  deviation  to  the  average  amplitude 

» ^ /dCaI  . (89) 

The  phase  fluctuations  were  measured  by  reeordiag  the  time  delay 
of  the  scattered  pulse  signal  relative  to  a continuous  wave  that  was 
displayed  on  the  asetlleseope . The  time  for  a wavelength  of  the  ev  signal 
was  measured  on  the  oscilloscope  (iQ.t  pace  for  95.8  kilsb  the  ratio  of 
these  two  times  multiplied  by  deg  gives  the  relative  pitas©  ^ in  degrees. 
The  average' phase  of  the  ten- different  measurements 
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was  used  to  calculate  the  variance 
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from  which  the  standard  deviation 

■ yffi*)  (98) 

was  obtained.  In  this  case,  the  standard  deviation  represents  the  phase 
fluctuations . 

The  reference  used  to  determine  the  scattering  coefficient  was 
obtained  experimentally  by  measuring  the  peak- to- peak  voltage  of  the 
pulsed  signal  reflected  from  a plane  pressure  release  surface  of  the  same 
dimensions,  32  in.  by  52  in.,  as  the  rough  surface.  Assuming  that  the 
pressure  of  the  underwater  sound  at  the  receiver  is  proportional  to  the 
voltage  output  f the  receiving  hydrophone,  one  defines  the  scattering 
coefficient  for  each  individual  measurement  as  the  ratio  of  the  scattered 
intensity  from  the  rough  surface  to  the  reflected  intensity  from  the  plane 
surface.  Thus,  the  experimental  scattering  coefficient  can  be  written  as 


where 

^ is  the  gracing  angle, 

\L  is  the  peak- to- peak  voltage  of  the  pulsed  signal  scattered  from 
the  rough  surface,  and 

V is  the  peak- to- peak  voltage  of  the  same  pulsed  signal  reflected 
P 

from  the  plane  surface. 
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The  average  scattering  coefficient,  which  is  the  correct  coefficient  to 
compare  with  the  theory  is 


«„<h 


(94) 


where 

i.  is  the  sample  number  (i  = 1,2, ...10),  and 
j is  the  grazing  angle  (j  a 2,4, ...80  deg). 


The  scattering  coefficient  expressed  in  decibels  is 

oe(dB)  = 10  logl0ce 


(95) 


CHAPTER  V 
RESULTS 


A.  Amplitude  Fluctuations 


The  experimental  results  for  the  amplitude  fluctuations  in  the 
specular  direction  for  a frequency  of  95.8  kHz  are  presented  in  Figs.  13 
and  14.  These  results  were  obtained  from  averages  of  readings  from  ten 
different  insonified  regions  of  the  rough  surface  by  using  Eqs.  (86) 
through  (89)  of  Chapter  IV.  The  theoretical  predictions  of  the  amplitude 
fluctuations  in  the  Fraunhofer  approximation,  given  by  Eq.  (24),  and  in 
the  Fresnel  approximation,  given  by  Eqs.  (70)  and  (72)  are  illustrated 
in  Fig.  13.  A bearawidth  of  9 deg  measured  at  the  half- power  points 


(-3  dB)  of  the  main  lobe  of  the  incident  beam,  as  shown  in  Fig.  8 of 
Chapter  IV,  was  used  to  calculate  the  dimensions  of  the  insonified  area. 
The  experimental  data  were  used  to  construct  the  dependence 

of  the  amplitude  fluctuations  on  the  Rayleigh  parameter,  2k  /\F”  ) sin  iff. 


Figure  13  pertains  to  measurements  made  on  the  model  surface  whose 
statistical  properties  were  characterized  by  a root-mean- square  relief, 

, of  0.091  inch.  From  the  graph  it  is  evident  that,  within 

fluctuations,  the  experimental  results  correspond  to  the  theoretical 


law  ’*«  the  Fresnel  approximation  for  Rayleigh  parameters  less  than  0.6. 


The  theoretical  law  in  the  Fraunhofer  approximation  of  zero  amplitude 
fluctuations  gives  incorrect  predictions  for  any  value  of  the  Rayleigh 
parameter.  The  experimental  amplitude  fluctuations  for  values  of  the 
Rayleigh  parameter  greater  than  one  could  not  possibly  be  explained  by 
means  of  the  theoretical  analysis,  as  the  region  of  applicability  of  the 
analytical  formulas  is  limited  by  the  condition  given  by  Sq.  (13)  that 


sin  f < 1. 


When  the  projector  and  receiver  were  aimed  at  each  other  over 


a distance  of  10  ft  the  amplitude  fluctuations  were  measured  to  be  less 


than  one-half  percent.  Thus,  changes  in  the  properties  of  the  water 


% 


FIGURE  13 

AMPLITUDE  FLUCTUATIONS 
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through  which  the  transmitted  and  received  pulsed  wave  propagated,  and 
inaccuracies  in  positioning  the  projector  and  receiver  contributed  less 
than  G.5$  to  the  amplitude  fluctuations  shown  in  Figs.  13  and  Ik. 

The  experimental  data  of  the  amplitude  fluctuations  depicted 
in  Fig.  14  are  the  outccme  of  combining  *;he  results  obtained  from  three 
different  model  surface,  I,  II,  and  III,  For  small  values  of  the  Rayleigh 
parameter  up  to  1.7  there  is  an  apparent  increase  of  the  amplitude 
fluctuations  from  10$  to  40$,  which  is  proportional  to  the  value  of  this 
parameter.  For  values  of  the  parameter  greater  than  2.5  the  growth  of 
the  amplitude  fluctuations  ceases  and  it  reaches  its  saturation  point. 

For  larger  values  of  the  parameter  there  is  no  further  increase  in  the 
fluctuations,  and  the  amplitude  fluctuations  disclose  a spread  of  about 
33  - 63$  in  the  values.  The  graph  shows  a well-defined  alternation  of 
the  amplitude  fluctuations  between  maximum  and  minimum  values.  The 


fluctuation  maxima  are  situated  at  2k  /(f  \ sin  a 2.5,  4.2,  5.6,  and  6.7 

K ' 

the  minima  at  2k  Air  / sin  i|*  0 3.4,  4.<),  and  6,4. 

The  measurements  illustrated  in  Fig.  14  verify  the  tentative 
results  obtained  by  Gulin  and  Malyshev  (1962)  who  scattered  sound  signals 
from  the  surface  of  the  sea.  Their  res' Jits,  which  were  stated  in  Chapter  I, 
show  that  the  amplitude  fluctuations  retch  a saturation  point  when  the 


Rayleigh  parameter  has  the  value  0.7  and  for  larger  Rayleigh  parameters 
the  fluctuations  experience  an  altemat  ion  between  maximum  and  minimum 
values  of  £5  to  ?0$. 

S.  Phase  fluctuations 


The  experimental  results  for  the  phase  fluctuations  in  the 
? locular  direct!.©®  for  a frequency  tzf  ^.S  kHa  are  presented  in  Figs.  13 
ani  16.  These  results  were  obtained  fi«m  averages  of  readings  from  ten 
different  tfisoanfied  regions  of  the  rough  surface'  by  using  B%«.  (§0) 
through  (9®)  of  Chapter  IV.  the  the  critical  predictions  of  the  prase 
fluctuations  in  the  Frawihofer  appratisatiea,  giver*  by  Sg.  (6S),  and  in 
the  Fresnel  approximation,  given  by  $$s.  (70#  and  (?b|,  are  depicted  in 
Fig.  l|.  $ 9 deg  besasfidih  was  used  i.a  the  calculation  for  <2^,  Gg,  add  I 
appearing  in  the  formulas . 
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through  which  the  transmitted  and  received  pulsed  wave  propagated,  and 
inaccuracies  in  positioning  the  Rector  and  receiver  contributed  less 
than  0.5$  to  the  amplitude  fluctuations  shown  in  Figs.  15  and  14. 

The  experimental  data  of  the  amplitude  fluctuations  depicted 
in  Fig.  14  are  the  outcome  of  combining  “he  results  obtained  from  three 
different  model  surface,  I,  II,  and  III.  For  small  values  of  the  Rayleigh 
parameter  up  to  1.7  there  is  an  apparent  increase  of  the  amplitude 
fluctuations  from  10$  to  40$,  which  is  proportional  to  the  value  of  this 
parameter.  For  values  of  the  parameter  greater  than  2.5  the  growth  of 
the  amplitude  fluctuations  ceases  and  it  reaches  its  saturation  point. 

For  larger  values  of  the  parameter  there  is  no  further  increase  in  the 
fluctuations,  and  the  amplitude  fluctuations  disclose  a spread  of  about 
55  - 65$  in  the  values.  The  graph  shows  a well-defined  alternation  of 
the  amplitude  fluctuations  between  maxim;  in  and  minimum  values.  Tne 

fluctuation  maxima  are  situated  at  2k  , 4f2\  sin  s 2.5,  4.2.  5.8,  and  6.7 ; 

/i\  , 7 , 

the  minima  at  2k  y \F^  ) sin  t a 5.4,  4.|),  and  6.4. 


The  measurements  illustrated  in  Fig.  14  verily  the  tentative 
results  obtained  by  Gulin  and  Malyshev  { 1962 ) who  scattered  sound  signals 
from  the  surface  of  the  sea.  Their  resets,  which  were  stated  in  Chapter  I, 
show  that  the  amplitude  fluctuations  reach  a saturation  point  when  the 
Rayleigh  parameter  has  the  value  0.7  aud  for  larger  Rayleigh  parameters 
the  fluctuations  experience  an  alternation  between  maximum  and  minimum 
valv.es  of  25  to  ?0$. 

B.  Fhas§  Pluet  u at i ana 

The  experimental  results  for  the  phase  fluctuations  in  the 
?t*?eular  direction  for  a frequency  of  *$,8  fcHs  are  presented  in  Figs.  13 
and  16.  These  results  were  obtained  ftom  averages  of  readings  free  ten 
different  insoar. fled  regions  of  the  rough  surface  by  using  113s.  (00) 
through  (98)  of  Chapter  IV.  The  theoretical  predict  lens  of  the  phase 
fluctuations  in  the  Fraunhofer  approx!  tat  ion,  given  by  (66),  and  in 
the  Fresnel  approximation,  given  by  ®?s.  (70)  and  are  depicted  in 
Fig.  13.  & 9 deg  beaswidth  was  used  !.n  the  calculation  for  £^,  and  I 
appearing  in  the  formulas . 


FIGURE  16 
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When  the  projector  and  receiver  were  aimed  at  each  other  over 
a distance  of  10  ft  the  phase  fluctuations  were  measured  to  be  less  than 
%,  This  small  fluctuation  in  phase  was  caused  by  the  inaccuracy  in 
moving  the  source  and  receiver  as  a rigid  body  relative  to  the  model 
surface,  and  also  was  brought  about  by  local  changes  in  the  properties  of 
the  water  through  which  the  acoustic  wave  propagated. 

The  dependence  of  the  phase  fluctuations  on  the  Rayleigh 
parameter  for  measurements  made  with  the  model  surface  that  has  a root-mean- 
square  relief  of  0.091  in.  is  shown  in  Fig.  15.  The  graph  shows  that  the 
theoretical  law  in  the  Fraunhofer  approximation  gives  predicted  values 
within  4$  of  the  values  obtained  using  the  Fresnel  approximation  for 
values  of  the  Rayleigh  parameter  less  than  one.  Over  this  range  of 
Rayleigh  parameters  it  is  evident  that,  within  15$  fluctuations,  the 
experimental  results  which  are  expressed  in  percent  relative  to  a radian, 
agree  with  both  theoretical  approaches.  The  phase  fluctuations  resulting 
from  the  irregularities  on  the  model  surface,  for  values  of  the  Rayleigh 
parameter  between  0.2  and  0.55  (theory  predicts  7 to  14$)  are  concealed 
possibly  by  the  fluctuations  caused  by  the  wave  traveling  directly  from 
the  source  to  the  receiver  and  the  wave  scattered  from  the  edge  of  the 
surface  interfering  with  the  wave  scattered  from  the  topograph  of  the  rough 
surface  (experimental  results  show  about  25$  relative  to  a radian).  With 
this  understanding  one  can  assert  that  both  theoretical  procedures  agree, 
within  6$  fluctuation,  with  the  experimental  results  for  values  of  the 
Rayleigh  parameter  less  than  0.6.  The  apparent  good  agreement  between  the 
experimental  results  and  the  theoretical  predictions  in  the  Fraunhofer 
approximation  for  values  of  the  Rayleigh  parameter  between  1.225  and  1.6? 
is  surprising  but  should  be  treated  as  purely  coincidental  since  the  theory 
is  not  valid  for  Rayleigh  parameters  greater  than  one. 

The  experimental  phase  fluctuations  illustrated  in  Fig.  l6  are 
the  consequence  of  incorporating  the  measurements  obtained  fro a three 
model  surfaces.  I,  II,  and  III.  The  dependence  of  the  phase  fluctuations 
on  the  Rayleigh  parameter  is  not  as  well  defined  as  the  amplitude 
fluctuations  on  this  parameter.  Nevertheless,  in  general  the  phase 
fluctuations  increase  with  the  Rayleigh  parameter  from  a low  value  of 


1 c.j>  to  a high  value  of  250$.  It  is  understood  that  the  percentage  for 
the  experimental  phase  fluctuations  is  percent  relative  to  a radian.  The 

fluctuation  maxima  are  situated  at  2k  sin  ^ » 1.5,  3.5,  4.6,  and  6.0 

the  minima  at  2k  sin  t « 1.8,  3.8,  5.3,  and  6.8.  Since  the  author 

is  unaware  of  any  other  experimental  investigation  of  the  phase  fluctuations 
of  an  underwater  acoustic  pulsed  wave  scattered  from  a randomly  rough 
surface,  no  comparison  can  be  made  and  these  results  should  be  considered 
tentative  and  in  need  of  further  verification. 

When  one  attempt,*-,  to  compare  the  experimental  phase  fluctuations 
with  the  corresponding  amplitude  fluctuations,  one  is  confronted  with  a 
problem.  The  issue  is  understood  when  it  is  realized  that  the  experimental 
amplitude  fluctuations  given  by  Eq.  (89)  represent  a relative  standard 
deviation,  whereas  the  experimental  phase  fluctuations  givrru  by  Eq.  ( 92) 
r** present  an  absolute  standard  deviation.  The.  experimental  amplitude 
fluctuations  in  Fig.  13  are  in  percent  relative  to  the  average  amplitude, 
whereas  the  phase  fluctuations  in  Fig,  15  are  in  percent  relative  to  a 
radian.  The  Fraunhofer  and  Fresnel  theory  for  the  amplitude  and  phase, 
fluctuations,  using  the  same  convention  for  the  ordinate  scales,  are 
presented  in  Figs.  13  and  15  in  percent. 

C.  Scattering  Coefficient 


the  experimental  results  for  the  scattering  coefficient  in  the 
specular  direction  for  a frequency  of  95,8  kK*  are  presented  in 
Figs.  17,  18,  19,  and  £0,  These  results  were  obtained  fyo*a  averages  of 
readings  from  ten  different  ineonified  regions  of  the  rough  yuriuce  by  - 
employing  (95)  through  (95)  of  Chapter  IV,  The  theoretical  predictions 
of  the  scattering  coefficient  in  the  Fraunhofer  approximation,  given  by 
Eq.  (79),  and  in  the  Fresnel  approximation,  given  by  M . (83),  are  alwm 
in  the  four  graphs.  The  dfasensions  of  the  iasonified  area,  and  i 

t>-*H . appear  In  the  foftaul&s  'mm  calculated  for  & 9 deg  beamridth 
1 ired  at  the  half- power ' potato  of  the  main  lobe  of  the  incident  beam. 
Figure  17  shows  t*  a, dependence  of  the  scattering  coefficient 
on  the  Bnyleigh  parameter  and  the  grating  angle,  .for  the  aodel  surface 

svrsss  belief  of  0,0$L  in,  and.- & correlation  length  of  t.55  iooheo. 


that  hue 
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Prom  the  graph  it  is  evident  that,  within  5 dB,  the  experimental  results 
correspond  to  the  theoretical  law  in  the  Fresnel  approximation  over  the 
complete  range  of  the  Rayleigh  parameters  from  values  of  0.20  to  1.85 . 

The  theoretical  law  in  the  Fraunhofer  approximation  gives  good  agreement 
with  the  experimental  results  up  to  a value  of  1.3  for  the  Rayleigh 
parameter.  But  for  larger  values  of  the  parameter  the  Fraunhofer  theory 
predicts  lower  values  for  the  scattering  coefficient  than  the  experimental 
results  show.  At  a Rayleigh  parameter  of  1.85  the  Fraunhofer  law  gives 
values  of  the  coefficient  that  are  4 dB  smaller  than  the  Fresnel  predictions . 

The  scattering  coefficient  measured  with  the  model  surface  that 
has  an  ms  relief  of  0.182  in.  and  a correlation  length  of  2.55  in.  is 
illustrated  in  Fig.  18.  Again,  it  is  evident  from  the  graph  that  the  theory 
in  the  Fresnel  approximation  predicts  the  experimental  values  better  than 
the  Fraunhofer  theory.  The  Fresnel  theory  predicts  the  experimental 
results  to  within  3 dB  for  Rayleigh  parameters  less  than  2.35.  For  larger 
values  of  the  Rayleigh  parameter  the  Fretmd  theory  predicts  values  within 
4.5  dB  of  the  measured  scattering  coefficient.  The  Fraunhofer  theory 
gives  good  agreement,  within  3 dB,  with  the  measured  values  of  the 
scattering  coefficient  for  values  of  the  Rayleigh  parameter  less  than 
1.25.  But  for  larger  values  of  this  parameter  the  theoretical  law  in 
the  Fraunhofer  approximation  deviates  notably  from  the  measured  coefficient, 
more  than  l6  dB  for  a Rayleigh  parameter  of  3.T. 


Figure  19  describes  the  measurements  and  predicted  values  of 
the  scattering  coefficient  for  the  model  surface  that  has  an  ms  relief 
of  0.164  ia.  and  a correlation  length  of  2.|J  laches.  As  la  the  two 
previous  eases  the  theory  derived  using  the  Fresnel., approximat ten  agree® 
with  the  experimental  result®  far  better  than  dees  the  Fraunhofer  theory. 
The  Fresnel  theory  1®  within  1$  dB  of  the  measured  values  of  the  scattering 
coefficient  whereas  the  Fraunhofer  theory  is  within  22  dB  for  tits  complete  . 
range  of  Rayleigh  parameters  shown  in  Fig.  If. 

A eoaparisoti  of  Fig®.  X?,  IS,  and  19  shswa  that  the  theoretical 
law  in  the  Fresnel  approximation  gives  progressively  worse  agreement  with 
the  expertee«t&l  restilts  a®  the  mo  relief  of  the  made!  surface  increases*,. 
The  reason  fre*  this  disparity  for  the  rougher  surfaces  ts  attributed . to 


the  breakdown  of  the  approximation,  ^ 


used  in  deriving  Eq.  (8) 


from  Eq.  (7)  in  Chapter  II.  As  shown  in  Table  I of  Chapter  IV  the  slope 
of  the  irregularities  on  the  surfaces  increases  at  the  same  rate  as  the 
rms  relief  of  the  first  three  model  surfaces.  Apparently  the  slopes  for 
model  surface  III,  results  depicted  in  Fig.  19,  are  too  great  to  justify 
the  replacement  of  differentiation  along  the  normal  with  differentiation 
along  the  negative  2- axis . 

Figure  20  illust; rates  the  experimental  measurements  and  the 
theoretical  predictions  of  the  scattering  coefficient  for  the  model 
surface  that  has  an  mis  relief  of  0.36k  in.  and  a correlation  length  of 
5.10  inches.  The  irregularities  on  this  model  surface  are  characterized 
by  the  same  slopes  as  the  model  surface  studied  in  Fig.  18.  The  theory 
in  the  Fresnel  approximation  predicts  values  for  the  coefficient  within 
7 dB  oi  the  measured  values  whereas  the  Fraunhofer  gives  values  within 


19  dB. 


FIGURE  20 
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CHAPTER  VI 
GENERAL  DISCUSSION 


A.  Scattering  Coefficient 

In  the  previous  chapter  the  dimensions  of  the  insonified 
area,  A^,  Ag,  and  £,  which  appear  in  the  theoretical  expressions,  are  cal- 
culated from  Eqs.  (B-l),  (B-2),  and  (B-3),  in  Appendix  C by  using  a 9 deg 
beamwidth  for  the  incident  radiation.  Also,  the  mathematical  analysis  was 
simplified  by  approximating  the  elliptical  scattering  region  with  a rectang- 
ular scattering  area  as  shown  in  Fig.  % of  Chapter  III.  The  9 deg  beamwidth 
was  chosen  because 

(l'  ..t  is  the  beamwidth  measured  at  the  half-power  points  (-3  dB) 

of  the  main  lobe  of  the  directivity  pattern  of  the  incident 
radiation  illustrated  in  Fig.  8 of  Chapter  IV,  a^d 
( '■?)  it  gives  the  best  agreement  of  thr;  theoretical  law  in  the 
Fresnel  approximation  with  the  experimental  scattering  coef- 
ficient measured  from  the  "smoothest''’  mode l surface,  rtns 
relief  of  0.091  in.,  as  sh^wn  in  Fig.  1. 

Figure  Si  illustrates  the  influence  the  -«ulue  of  the  beamwidth, 
and  htnct  the  dimensions  that  the  insonified  area,  ua;  on  the  theoretical 
predict lens.  At  the  Rayleigh  parameter  of  1.8  tfefe  8 dig  beamwidth  results 
in  predicted  values  1.8  dB  smaller  than  the  Qbsewed  scattering  coefficient, 
the  f deg  beamwidth  predicts  0.?  dl  larger  val>uef  > the  10  deg  beamwidth 
predicts  k.6  di  larg.tr  valuta,  and  the  12  deg  beamwidth  gives  values  that 
art  J.t  dl  larger  than  the  measured  value®  of  the  sc&ttiring  coefficient. 

fhe  effect  that  a change  in  the  beamw.Lath  of  the  incident  radiation 
hag  on  tht  Fresnel  theory  for  the  roughest  surface,  rms  relief  of 

0.36&  in.,  is  illustrated  in  Fig.  riZ.  With  reference  to  a value  of  ?,1  f@r 
the  Rayleigh  parameter  it  is  evident.  m«  the  graph  that  a $ deg  beamwidth 
results  in  a the@retie.ai  value  for  1 he  scattering  eoeffieient  that  it  i£  dl. 
smaller  than  the  measured  coefficient,  an  8 deg  beamwidth  predicts  a 13. 5 dB 
smaller  value,  a 9 deg  besawidth  gives  a $ 'dB  smaller  value.  10  deg  besa- 
vtdth  gives  a t.5  dB  lower  >alue,  mi  a 12  deg  beaswldth  results  in  a ^ dB 
lower  value'. 
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Even  though  the  12  deg  hem  results  In  predicted  values  that  give 
best  agreement  with  the  observed  scattering  coefficient  for  the  roughest 
surface,  a 9 deg  beamwidth  was  chosen  since  this  selection  agrees  best  with 
the  experimental  results  of  the  smoothest  surface  shown  in  Fig.  21. 

Table  I in  Chapter  IV  shows  that  the' slopes  of  the  irregularities  of 
model  surface  III  are  four  times  as.  great  as  the  slopes  on  model  surface  I, 
the  smoothest  surface.  One  would  expect  the  theory,  which  is  subject  to 
the  approximation  <D/du  ~ -o/c)z,  to  give  its  best  predictions  for  the  experi- 
mental results  obtained  with  the  model  surface  that  is  characterized  by  the 
least  slopes  of  the  irregularities-. 

Figure  23  is  a plot  of  Kqs.  (53)  and  (5*0  of  Chapter  II, 
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which  is  proportional  to  the  intensity  in  the  Fresnel  approximation  of  a wave 
scattered  fire®  a plane  surface  (the  denominator  of  the  scattering  coefficient). 
The  change  in  the  intensity  with  change  in  the  beamwidth  is  the  result  of 

* 

the  prtsenee  of  the  Fresnel  integrals,  C and  §»  In  the  expression  for  p^p  . 
This  relation  has  been  discuss  ad  thoroughly  by  Leittr  ClfH).  A result  ©£ 
his  study  was  that  tht  tsr®  K tends  t©  unity  in  a fluctuating  manner  as- 
the  number  of  Fresnel  senes  enclosed  within  the  rectangular  Insenlfied  area 
increases*  Figure  f|  exhibits  tfe.g  behavior  very  clearly  for  a®  the  beam- 
width  ©f  the  incident  radiation  r increased  without  changing  the  source 
distance  (8^),  receiver  distance  (R^),  and  the  frequency,  the  number  ©f 

Fresnel  ecnes  enclosed  in  the  insenifiei  area  increases.  The  graphs  shew 

■ - - ~ «... 

that  as  the  beamwidth  increases  frees  6 deg  t®  It  deg  the  term  p^p  apf-reaches 
1/(^0  t which  is  tht  sts*; tight  line  with  the  value  X 10**^. 

Figure  I»b  in  Appendix  1 shew©  fiat  a p deg  beamwidth  produces  .an  insenified 
area  that  encloses  only  one  Fresael  sene,  k 6 dsg  beamwidth  will  inss&ify  ' 
an  area  less  than  ®ne  Fresnel  sene  whereas  a IS  deg  beamwidth  will  inssnify 


approximately  the  first  two  Fresnel  zones..  Thus,  the  solution  for  the 
intensity  reflected  from  a plane  surface  calculated  in  the  Fresnel  approx-, 
imation  for  an  active  scattering  region  that  is  rectangular  will  become 
steady  as  the  value  K = 1 is  approached,  and  in  turn  the  scattering  coef- 
ficient will  become  stable. 


Leizer  (1966)  has  shown  also  that  the  stability  of  the' intensity 
reflected  from  a plane  surface  is  strongly  governed  by  the  shape  of  the 
insonified  area  used  in  the  calculations  in  the  Fresnel  approximation.  His 
study  has  shown  that  when  an  insonified  area  of  elliptical  shape  is  used 

•ifc 

in  the  calculation  of  p p the  term  K will  fluctuate  between  constant 

o o 

limits  0 and  it.  The  indicated  minima  (K  = 0)  and  maxima  (K  = L)  we  obtained 
when  an  even  or  odd  number  of  Fresnel  zones,  respectively,  exactly  fit  into 
the  elliptical  region.  Thus,  if  an  elliptical  insonified  area  had  been  used 
in  the  derivation  of  the  scattering  coefficient  in  this  paper,  an  infinite 
value  for  this  coefficient  would  have  resulted  whenever  a beamvid^h  was  chosen 
that  insonified  an  even  number  'f  Fresnel  senes.  The  infemation  of  primary 
importance  obtained  from  Leizer’ s study  was  that  the  solution  for  the 
reflected  intensity  would  not  stabilise  as  the  number  of  Fresnel  gene*  encom- 


passed by  the  elliptical  insonified  aim  increased,  whereas  the  solution 
would  stabilise  as  the  number  of  Fresnel ’sen^t  enclosed  by  the  rectangular 
insonified  area  increased. 


Tills  phenomenon  dees  net  ©oeu?  experimentally  because  the  insonified 
area  is  not  wsll-definti  as  Fig.  §•  shows.  Consequently,  tht  fluctuating 
behavior  @f  the  solution  in  the  Fresnel  approximation  ceuld  be  ©iiednsted  if 
better  approximations  of  the  transducer  beam  patterns  ©ou\d  be. Inserted''  . 
into  the  theory.  Merten  and  Huir.  ilffi?)  have  used  the  expression 


to*  daSorifei  the  effect  ®f  the  rs&intieft  pattern  @f  the  source  on  the  scat- 
tered intensity,  ffeia  relation  describes  a transducer  pattern  whose  mjftt 
lob®  inssnifies  an  nrer  of  • elliptical  shap*  on  the  suri^- , 45*4  that 
does  net  have  any  <sinar  leVse.  Their  ealculatiene  were  &s4e  .1®.  the  Fraunfesfer 
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approximation.  It  is  conceivable  that  the* theory  presented  in  the  current 
paper  ccxLd  be  improved  if  the  above  approximation  of  the  radiation  pattern 
of  the  source  could  be  inserted  into  the  theoretical  development  in  the 
Fresnel  approximation.  However,  one  may  encounter  insurmountable  analytical 
difficulties. 

The  scattered  intensity  calculated  by  using  the  statistical 
parameters  of  model  surface  I and  model  surface  III  is  plotted  versus  the 
grazing  angle  in  Figs.  24  and  25,  respectively.  The  term  pp*  is  the  numer- 
ator of  the  scattering  coefficient  and  represents  the  intensity  except  for 
a factor  of  l/2Z  where  Z is  the  specific  acoustic  impedance  of  the  medium.  - 
The  graphs  show  t-  there  is  only  a small  change  in  the . intensity  as  the 
beamwidth  varies  ij.om  8 deg  to  12  deg  compared  with  the  fluctuating  behavior 
of  the  formulas  derived  using  the  plane  surface  as  depicted  in  Fig.  2?. 

Thus,  the  variability  of  the  scattering  coefficient  with  beamwidth,  as 
illustrated  in  Figs.  21  and  22,  is  attributed  mainly  tc  the  fluctuating 
behavior  of  the  theoretical  expression  for  th?  intensity  reflected  from  a 
plane  surface,  the  denominator  of  the  scattering  coefficient. 

B.  Fraunhofer  and  Fresnel  Approximations  for  the  Intensity 

The  theoretical  scattering  coefficient  in  the  Fraunhofer 
approximation  gave  progressively  worse  predictions  for  the  observed  scat- 


tering  coefficient  as  the  rme  relief  y and/or  as  the  Rayleigh  parameter 

of  the  rough  surface  increased.  This  deviation  between  the  Fraunhofer  theory 
and  the  experimental  results  was  evident  in  Figs.  17,  3 8,  and.  19  of  Chapter  V. 
Thus,  the  Fraunhofer  approximation  predicts  incorrect  values  for  the  scat- 
tering coefficient  as  defined  by  Eq.  (48)  of  Chapter  II,  except  possibly  for 
the  very  "smoothest"  surfaces. 

This  phenomenon  is  reversed  when  only  the  intensity  of  the  wave 
is  considered.  From  Eq.  48  the  intensity  is 


FIGURE  24 

RELATIVE  INTENSITY 
HE  SPECULAR  DIRECTION 
IE  RELEASE  RANDOM  SURFACE 


FIGURE  25 

RELATIVE  INTENSITY 
IN  THE  SPECULAR  DIRECTION 
PRESSURE  RELEASE  RANDOM  SURFACE 


In  the  Fraunhofer  approximation  poPQ  is  given  by  Eq.  (C-5)  in  Appendix  C, 
and  a is  given  by  Eq.  (79)  in  Chapter  III.  These  relations  give  for  the 
Fraunhofer  intensity 


7b 


, 2. 2, . A ,2  . 2 

/A  ^ + A2)  Sin  * _-g 

\"v  FR  " ? £ * e 


2pvit  R2q 


1 + 


2 

a it 


m=oo 


V 


m 


T(Z  + V 2 

I c , m.m 
m=l 


(97) 


In  the  Fresnel  approximation  p p is  given  by  Eq.  (C-12)  in  Appendix  C, 

0 0" 

and  a by  Eq.  (83)  in  Chapter  III.  These  egressions  give  for  the  Fresnel 
s 

intensity 
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The  ratio  of  the  Fraunhofer  to  the  Fresnel  intensity  is 
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Figure  ?6  is  a plot  of  this  ratio  in  decibels  given  by  the  relation 


dB  = 10  log. 


for  the  plane  surface  and  the  four  model  surfaces  listed  in  Table  I of 
Chapter  IV.  A 9 deg  beamwidth  -was  used  in  the  calculations.  At  a grazing 
angle  of  80  deg  the  ratio  for  the  plane  surface  is  12  dB,  for  surface  I 
the  ratio  is  8 dB,  for  surface  II  the  ratio  is  0.7  dB,  and  for  sv.rfaceg 
III  and  IV  the  ratio  is  less  than  0.2  dB.  From  the  graph  it  is  evident  that 
for  grazing  angles  greater  than  30  deg  the  intensity  calculated  in  the 
Fraunhofer  approximation  is  within  1 dB  of  the  intensity  calculated  in  the 


FIGURE  26 

RATIO  OF  FRAUNHOFER  TO  FRESNEL  INTENSITY 


Fresnel  approximation  for  the  roughest  model  surface.  But  for  the  plane 
surface  the  two  theories  differ  "by  as  much  as  13  dB  for  grazing  angles 
greater  than  J>0  deg.  Also,  the  effect  that  the  correlation  length  has  on 
the  approximations  is  revealed.  As  the  correlation  length  of  a rough  surface 
decreases,  the  intensity  derived  using  the  Fraunhofer  approximation 
approaches  the  values  predicted  by  the  Fresnel  approximation.  Therefore, 
the  theoretical  law  for  the  intensity  scattered  from  a rough  surface  cal- 
culated in  the  Fraunhofer  approximation  will  give  poor  agreement  with 
the  observed  intensity  for  the  smoother  surfaces  and  good  agreement  for  the 
rougher  surfaces. 

Figure  27  is  a plot  of  the  same  ratio  vs  the  Rayleigh  parameter 
for  model  surfaces  I,  II,  and  III.  For  Rayleigh  parameters  greater  than  ^ 
for  the  numerical  example  used  throughout  this  paper,  the  Fraunhofer 
approximation  will  give  values  for  the  intensity  that  are  within  1 dB  of 
the  values  predicted  by  the  theory  using  the  Fresnel  approximation.  Thus, 
the  choice  to  use  either  the  Fraunhofer  or  Fresnel  approximation  to 
study  the  intensity  scattered  from  a surface  will  depend  upon  the  roughness. 
A study  of  the  scattering  coefficient  will  require  the  Fresnel  approxi- 
mation in  every  case  except  for  a very  smooth  surface  where  the  Rayleigh 
parameter  is  less  than  one. 

C . Application  of  Model  Studies  to  Ocean  and  Shallow  Water  Studies 

The  slant  line  distance  from  the  source  to  the  scattering 
surface  in  the  model  studies  was  5 ft.  A common  distance  from  source  to 
the  sea  bottom  for  ocean  studies  is  5000  ft,  which  requires  a factor  of 
1000  to  scale  the  model  studies  to  an  actual  situation  in  the  ocean. 

Shallow  water  studies  are  conducted  in  many  instances  at  a depth  of  100  ft. 
In  this  case  a scale  factor  of  20  is  needed.  The  type  of  source,  diameter 
of  source,  and  frequency  of  operation  needed  for  these  studies  to  be 
directly  applicable  are  listed  in  Table  II.  Also,  the  rms  relief  and  the 
correlation  length  of  the  model  surfaces,  listed  in  Table  I of  Chapter  IV, 
should  be  scaled  to  conform  with  the  shallow  water  and  ocean  studies. 


CHAPTER  VII 
CONCLUSIONS 

The  theoretical  approach  of  Gulin  (1962)  has  been  used 
successfully  to  predict  the  observed  amplitude  and  phase  fluctuations 
of  an  acoustic  wave  scattered  in  the  specular  direction  from  a rough 
surface  that  is  described  statistically  by  an  exponential  covariance 
function.  The  theoretical  expressions  were  limited  to  values  of  the 
Rayleigh  parameter  less  than  one.  The  comparison  between  experimental 
and  theoretical  results  showed  the  necessity  of  the  Fresnel  theory  to 
describe  the  amplitude  fluctuations.  Fraunhofer  and  Fresnel  approximations 
gave  equally  good  agreement  with  the  experimental  phase  fluctuations. 

The  experimental  amplitude  and  phase  fluctuations  were  presented 
for  values  of  the  Rayleigh  parameter  less  than  7*^*  The  observed  amplitude 
fluctuations  in  the  specular  direction  for  the  model  studies  showed  the 
same  behavior  as  the  measurements  made  xa.  the  sea  by  Gulin  and  Malyshev 
(1962),  who  limited  their  study  to  values  of  the  Rayleigh  parameter  less 
than  3.5.  The  author  is  unaware  of  any  other  messurememts  of  the  phase 
fluctuations  of  underwater  acoustic  waves  scattered  from  a rough  surface. 
Since  accurate  measurements  of  the  phase  fluctuations  are  difficult, 
further  experimental  studies  are  desirable. 

The  theoretical  development  hud  to  be  extended  to  the 
Fresnel  approximation  to  obtain  values  for  the  scattering  coefficient 
that  were  in  good  agreement  with  the  measured  values.  Even  though  the 
theoretica.1  model  produced  values,  which  fluctuated  notably  whenever 
only  a few  Fresnel  zones  were  insonified,  the  theory  in  the  Fresnel 
approximation  gave  significantly  better  agreement  with  the  experimental 
scattering  coefficient  than  the  Fraunhofer  approximation  for  the  four 
randomly  rough  model  surfaces  used  in  this  work. 

The  theory  of  rough  surface  scattering  presented  in  this 
paper  used  the  following  assumptions  or  simplifying  procedures: 

1.  It  was  assumed  that  all  the  radiation  that  impinged  on  the 
rough  interface  was  either  reflected  or  scattered. 

2.  Multiple  scattering  was  neglected. 
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3.  Shadowing  effects  were  neglected. 

4.  The  radius  of  curvature  of  the  scattering  elements  was  taken 
to  be  much  greater  than  the  wavelength  of  the  incident 
radiation.  This  restriction  along  with  assumptions  1,  2, 
and  3 allowed  the  field  at  each  point  of  the  surface  to 

be  represented  as  the  sum  of  the  incident  wave  and  a wave 
reflected  from  the  plaue  tangent  to  the  surface  at  the 
given  point.  These  statements  constitute  the  Kirchhoff 
approximation. 

5.  The  distances  of  the  source  and  the  receiver  from  the 
scattering  surface  were  assumed  to  be  much  greater  than  a 
wavelength  of  the  incident  radiation. 

6.  The  irregularities  were  assumed  to  have  sufficiently  small 
slopes  so  that  differentiation  along  the  normal  to  the 
surface  can  be  approximately  replaced  by  differentiation 
along  the  z-axis.  Beckmnn  (19&3)  avoided  waking  this 
assumption  by  inserting  into  his  theory  the  complete  ''local'* 
scattering  geometry.  However,  the  additional  terns  obtained 
by  making  this  insertion  reduced  to  uni 1 y whenever  the 
problem  was  restricted  to  the  special  ease  of  specular 
scattering. 

?.  Edge  effects  caused  by  the  finite  site  of  the  model  surfaces 
were  neglected  in  the  theory. 

8.  The  amplitude  of  excitation  of  the  incident  radiation  was 
assumed  to  be  practically  constant  over  the  active 
scattering  region. 

9.  The  theory  for  the  amplitude  and  pirns©  fluctuations  was 
restricted  to  values  of  the  Rayleigh  parameter  less  than  one. 

10.  The  probability-density  function  that  governs  the 
simultaneous  occurrence  of  the  elevations  on  the  rough 
surface  was  assumed  to  be  a bivariate  O&ussian  distribution. 

11.  The  rough  surface  was  assumed  to  be  randomly  rough, 
isotropic,  and  stationary  in  the  wide  sense. 

12.  All  calculations  were  restricted  to  the  specular  direction. 


13.  The  Fresnel  approximation  of  neglecting  cubic  and  higher 
ordered  terms  was  employed. 

14.  The  active  scattering  area  was  assumed  to  be  large  enough 
so  that  many  surface  irregularities  were  enclosed  within 
this  region. 

13.  In  the  application  of  the  theoretical  formulas  to  an 

actual  situation  it  was  assumed  that  the  active  scattering 
region  was  the  insonified  area  irrespective  of  the  number 
of  Fresnel  zones  insonified. 

l6.  The  beam  pattern  of  the  energy  radiated  from  the  source 
was  approximated  by  a cone  that  had  a beamwidth  measured 
at  the  half- power  points  (-3  dB)  of  the  actual  transducer 
radiation  pattern. 

IT.  The  elliptical  insonified  area  resulting  fro®  the 

intersection  of  the  conical  beam  of  incident  radiation 
with  the  scattering  surface  was  approbated  by  a rectan- 
gular active  scattering  region. 

Although  this  list  of  assumptions  is  lengthy,  they  are  utilised 
by  nearly  all  workers  in  this  field.  The  substantiated  agreement  between 
experiment  and  theory  shows  the  assumptions  are  physically  acceptable. 

The  major  innovation  of  the  present  study  is  assumption  13.  Although 
the  Fresnel  assumption  has  been  utilised  in  a few  theoretical  studies, 
nowhere  lias  a systematic  study  been  made  of  the  importance  of  the  Fresnel 
terms. 

in  future  studies  it  may  be  possible  to  withdraw  seme  of  these 
assumptions  at  tlie  cost  of  mathematical  simplicity.  For  example,  it  may 
be  possible  to  remove  assumption  l€  by  using  a better  approximation  for 
the  radiation  pattern  of  the  source  such  as  was  discussed  by  Horton  and 
Muir  (19&?).  It  Is  very  lirely  that  restriction  12  eorid  be  generalised 
somewhat  to  include  the  ease  of  any  receiving  angle  independent  of  the 
g rasing  angle  in  the  plane  of  incidence. 


the  mmmwm  md  eresnel  apphqxb^tions 


Figure  3 depicts  the  coordinate  system  with  the  origin  at  0. 
The  geometry  of  the  scattering  of  acoustic  waves  from  a randomly  rough 
surface  with  the  origin  at  O'  is  illustrated  in  Fig.  A-l.  As  in  Fig.  3 
the  source  Q and  the  receiver  A are  confined  to  the  x-z  plane.  This 
translation  of  the  origin  along  the  x-axis  is  expressed  as 
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and  for  the  specular  direction  can  he  written  as 


x . = x - 

1 «!  * 


(A-l) 


With  reference  to  Fig.  A-l  the  distance  can  be  expanded  in 
terms  of  R^.  Wttn  the  angle  between  r and  denoted  as  the  law 
of  cosines  gives 
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This  equation  eta  be  expressed  in  terns  of  the  direeU@a&i  angles • 
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In  a abailar  way  R can  be  written  in  terns  if  to  give 
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Adding  Eijs.  (A-3)  and  (A-6),  one  obtains 
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higher  order  tors? 


In  the  integrals  appearing  in  Chapter  IS  the  8 dependency  is 
taken  ml  ©f  as  illustrated  in  &js.  (U)  an*  (15)  e©  that  one  is 

interested  in  ^ * BJ.  Thus,  for  * * 0,  ^ * rj  and  8g  * and  g*.  (A-7) 


reduces  to 


.r_-.  f'ixawwt* *^«^?***.w* 
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g/ cos  lir  COS  * \ 

+ x.yn s rj-  + quart  ic  and  higher  order  terms 

l 2R2B  2R10  i 


The  Fraunhofer  approximation  is  the  neglect  of  the  quadratic  and 
higher  order  texmo  in  Eq.  (A-8).  The  quadratic  terms  are  insignificant 
whenever 


. t/x  , j_ 

P ^10  W W>  8Ra>, 


(A-9) 


where 


k ia  the  wavenumber, 


&*.  is  the  max iaum  dimension  of  the  active  scattering  reelun  in  the 
♦x^-direetion,  and 

I is  the  maximum  dimension  of  the  active  scattering  region  in  the 
*y— direction  and  t are  discussed  in  detail  in  Appendix  B). 

it  is  seen  from  solving  Kq . (A-9)  fur  that  the  receiver  will  be  in 
the  Fraunhofer  cone  (farfield)  when 


k sift  ifg  * t*J 
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The  Fresnel  approximation  is  the  neglect  of  the  cubic  and 
higher  order  terms  in  Eq.  (A-8).  The  receiver  uill  be  situated  in  the 
Fresnel  /’.one  when 
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Solving  for  R^,  one  obtains 
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Thus,  the  receiver  lies  within  the  Fresnel  zone  when 
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The  prnurth«fsr  approx imati-jn  t«  fjq.  ran  be  written  ft 


8i  • ^ R10  * "so 


ex. 


{A-lhi 


ff,  “ es» 

1 t 


e - cos 


( A- 1 > ) 


M^iiwrfMitiitfi'iiiiri'  iiif'ifvriiitftiii 


The  Fresnel  approximation  to  Eq.  (A-8)  can  be  written  as 


where 


^ + R2  = R10  + R20  + CXi  + Rg  X1  + R y 


R gRioR?o 

R10  3X11  ^2  4 R20  S^q 


, and 
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For  the  specular  direction 


\ s * .» 


and  Eq.  (A-8)  reduces  to 


*1  4 R2  “ ^0  * ^20  4 R ( *1  sin  ^ 4 y ] 


+ quortic  and  higher  order  terns 


Rote  that  the  linear  and  cubic  teres  V v<x  vanished  so  that  the  relation 
la  exact  through  cubic  terns  Tor  the  travel  uppraximation  In  the  specular 
direction.  The  Fraunhofer  approximation  in  the  specular  direction  is 
staple 

^ 4 rv2  9 R10  4 * (*-»>) 

When  the  Fraunhofer  approximation,  Eq.  (A-SO)  er  the  Fresnel 
approximation.  Es*.  (A-19),  is  made  in  Eqs.  (14),  (IT),  (§0),  or  (tl)  of 
Chapter  il  only  the  approximations 


®g  * 
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The  Fresnel  approximation  to  Eq.  (A-8)  can  be  written  as 


where 


b;  * r2  ■ \0  * “2o  * ch  * r 4 * I j2  ■ 
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?R10R?0 


R10  Sin2  ^2  + R20  si*  h 
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2R10R20 


R10  + 
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(A-18) 


For  the  specular  direction 


^ = 0?  = 'll 


and  Eq.  (A-8)  reduces  to 


R1  + H2  * **10  * R2Q  + | (*1  3i^  ''  + * 

+ quurtie  and  higher  order  terns 


(A-19) 


Note  that  the  linear  and  cubic  '<^se  have  vanished  so  that  the  relation 
La  exact  through  cubic  terns  f or  tn**  Fresnel  approximation  in  the  specular 
direst  ion.  The  Fraunhofer  approx i&at ion  in  the  specular  direction  is 
staple 
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When  the  Fraunhofer  approx l sat  ion,  gq. 
approx iaatiofi.  Eq,  (A- 19),  i»  s&de  in  Bqs.  (1W)# 
Chapter  tl  oai>‘  the  approximations 


and 


0)  or  the  Fresnel 
, (§0),  or  (a)  of 


«i  = ho 


9 


(iVSl) 
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are  used  tn  all  the  terns  of  the  Integrand  except  in  the  tern  k(R^  + R^) 
that  represents  the  phase.  What  is  implied  by  these  approximations  is  that 
at  distances  from  the  uneven  surface,,  large  compared  to  the  dimensions 
of  the  active  scattering  region,  the  amplitude  of  the  pressure  wave 
produced  by  any  particular  surface  element  differs  only  slightly  from 
that  produced  by  any  other  element.  On  the  other  hand  the  relative  piiase 
of  the  pressures  produced  at  the  receiver  by  any  two  surface  elements 
depends  on  the  differences  in  the  distance  of  the  two  elements  from  the 
receiver,  and  for  distant  points  this  difference  is  practically  inde- 
pendent of  P.;  and  R* . As  a consequence,  one  must  use  tie  entire 

O.  C. 

expressions  gion  in  Eqs.  (A- 19)  and  (A-20)  in  the  term  k(R|  + R^) 
appearing  in  the  integrands  of  Eqs.  (14),  (17),  (20),  and  (21). 

Consider  a numerical  calculation  of  the  Fraunhofer  and  Fresnel 
distances,  R^q,  as  given  by  Eqs.  (A- 10)  and  (A- 12).  The  source  is  placed 
60.0  in.  (R^)  from  the  surface.  With  the  source  and  :*eceiver  maintained 
in  the  specular  direction,  the  grazing  angle  of  the  source  is  varied. 

It  is  desired  to  know  where  the  receiver  must  be  placed  relative  to  the 
surface  (R^)  so  that  the  experimental  data  can  be  compared  with  theory 
that  was  developed  using  either  the  Fraunhofer  or  the  Fresnel  approximation. 
These  positions  of  the  receiver  are  referred  to  as  the  Fraunhofer  field 
sons  or  the  Fresnel  field  so.no,  whichever  the  case  may  be.  The  calculations 
are  made  for  a freqeuncy  of  9KM  kHz  and  for  a 9 deg  beaawldth  of  the 
source.  The  active  region  of  the  surface  Is  taken  to  be  the  inaonified 
area  whose  dimensions  are  calculated  from  Sqs.  (B-l),  (3-2),  and  (8-5)  in 
Appendix  8.  The  results  for  this  plane  surface  example  are  illitstrated 
in  Fig.  A- 2,  The  calculations  show  that  it  is  physically  Impossible  to 
place  the  receiver  in  the  Fraunhofer  field  cone  for  specular  reflection 
free  a Plane  surface.  The  receiver  aunt  be  placed  at  distances  such 
greater  than  44  in.  (tL,)  from  the  pressure  release  plane  surface  at 
§£ J kits  in  order  that  it  resain  in  the  Fresnel  field  gone  »t  all  receiving 
a>igles . The  abrupt  change  in  the  boundary  of  the  Fresnel  field  gone  at 
22  deg  receiving  angle  is  the  result  of  the  finite  dimensions  of  the 
reflecting  surface,  that  is,  is  held  constant  * 16  in.  for  gracing 
angles  less  than  22  deg. 
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In  Chapter  VI  it  is  shown  {**«  Fig.  27)  that  as  the  roughness 
(the  Rayleigh  jvxrameter  which  for  a plane  surface  Is  sero)  iacresaMfS, 
the  Fraunhofer  field  zone 'coves  in  from  infinity  to  see*  finite  distance 
from  the  rough  surface.  Figure  27  of  Chapter  VI  shows  for  the  numerical 
example  cited  tiiat  for  values  of  tne  Rayleigh  parameter  greater  than  A 
the  Fraunhofer  approximation  predicts  quantities  that  are  well  within 
1 dB  of  the  values  obtained  using  the  Fresnel  approximation.  The 
interpretation  of  this  behavior  is  that  for  very  rough  surfaces  the 
irregularities  on  the  surface  act  as  individual  scatterers.  Thus,  one 
should  not  insert  into  Eqs.  (A- 10)  and  (A-22)  the  dimensions  of  the 
active  scattering  region,  6^  and  f,  but  should  substitute  the  dimensions 
of  the  surface  Irregularities.  For  example,  if  one  teak  the  correlation 
length  a (2.55  In.  for  three  of  the  model  surfaces)  for  the  maximum  x and 
y dimensions,  Eq.  (A- 10)  would  give  for  the  previously  cited  example  for 
a grazing  angle  of  AO  deg 

Rjjq  » 15  m. 


for  the  Fraunhofer  field  zone. 

In  summary,  for  Rayleigh  parameters  much  less  than  one  the 
dimensions  of  the  active  scattering  region  are  used  in  Eqs,  (A-1Q)  and 
(A-lt).  When  only  a few  Fresnel  tones  on  the  surface  are  Immnlflecl  the 
dimensions  of  the  active  scattering  region  are  taken  to  be  the  dimensions 
of  the  insert i fled  area  (see  Appendix  B).  For  Rayleigh  parameters  much 
greater  than  one  the  dimensions  of  the  irregularities  on  the  surface  are 
used.  A satisfactory  approximation  of  the  sUe  of  the  irregularities  is 
the  correlation  length,  a. 


ACTIVE  SCATTERING  REGION 


What  regions  of  the  scattering  surface  are  the  most  important  in 
contributing  to  the  total  field  at  a given  receiving  point  when  the  surface 
is  illuminated  by  a source  at  a given  point?  Beckmann  (196?)  has  stated 
that  the  answer  is  "the  first  few  Fresnel  zones"  located  on  the  scattering  v 
surface.  For  an  omnidirectional  source  one  would  use  the  dimensions  of 
the  first  few  Fresnel  zones  as  the  limits  on  the  integrals  appearing  in 
Chapter  II . 

In  the  experimental  work  a directional  source  was  used,  which  has 
a directivity  pattern  illustrated  in  Fig.  B-l.  Only  a finite  portion  of  the 
scattering  surface  is  insonified.  The  insonifiod  area  will  be  defined  as 
the  intersection  of  the  half-power  points  (-3  dB)  of  the  main  lobe  of  the 
directivity  pattern  with  the  scattering  surface.  It  is  of  interest  to 
determine  both  the  dimensions  of  tbs  first  few  Fresnel  zones  and  the  dimen- 
sions of  the  irradiated  surface. 

The  length  of  the  irradiated  surface  is  determined  by  the  angular 
width  of  the  directivity  characteristic  of  the  transmitter,  the  slope  of  its 
axis,  tu.d  the  distance  from  the  transmitter  to  the  scattering  surface. 

The  dimensions  of  the  insonified  area  can  be  derived  in  the  following  manner. 
From  Fig.  B-i 

A » l80  deg  »■  % - V, 

where  ifn  is  the  gracing  angle  and  9 is  the  beaawidth  measured  at  the  half- 
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(B-2) 


With  the  projector  confined  to  the  x^-z  plane  the  linear  dimensions  of  the 
insonified  area  in  the  y direction  is  calculated  easily  from  the  right 
triangle  in  Fig.  P-2a  to  give 

£ = R1Qtan|  (B-3) 

The  insonified  area,  shown  in  Fig.  B-2b,  is  an  ellipse  with  a 
semimajor  axis  given  by 

“ * i ‘h  + V - 

and  a semiminor  axis  given  by 

l(\  + &,) 

b = — -~1— - ===== 

+ £^\ 

The  area  of  an  ellipse  is. nab,  so  that 

rt/(A.  * A-,)2 

area  = — -j'  — — =: ■ (B~6) 
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Kerr  ( 1951)  gives  formulas  for  the  calculation  of  the  dimensions 
of  the  Fresnel  zones.  The  geometrical  arrangement  of  source  and  receiver 
for  these  calculations  is  .shown  in  Fig.  B-3.  The  origin  lies  on  the 
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reflecting  plane  directly  below  the  source.  The  receiver  is  in  the  specular 
direction.  The  locus  of  all  points  in  the  x - y plane,  from  which  the  sec- 
ondary radiation  arrives  at  the  receiver  with  a constant  phase  difference  6 
with  respect  to  the  direct  radiation  along  r,  is  given  by  the  relation 

R[  + R2  * r = & 

■"he  path  difference  6 via  the  geometrical  point  of  reflection  is 

R10  + R20  " r = &o 

The  nth  ellipse  (Fresnel  zone)  is  determined  from 


6 - 6 = n % 
n o 2 


or 

6n  = *10 +' R20  - r * " ! (B'T) 

where  \ is  the  wavelength  of  the  incident  radiation.  Expressing  the 
ellipses  in  Urms  of  L,  6^,-z^,  &hd  006  ottains  for  the  x coordinate 
of  the  center  of  the  ellipse 


where  the  origin  w-jpicted  in  Fig.  B-|  has  been  transformed  to  the  point  of 
geometrical  reflection  by  means  of  Eq,  (A-i)  in  Appendix  A.  The  semiminer 
axis  is 


(e-9) 
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■anil  the  semimajor  axis  is 


a - b 
n n 
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The?  x intercepts  relative  to  tM 'geometrical  point  of  reflection  are 

A.  = a - x . and 

in  n on 

- a + x 
chi  n on 

?rom  Eq.  (B-5)  the  y intercepts  are 


*n  " +- 2T  /1(V 2 " 3 4 + 4 


(B-ll) 

(B-12) 

(B-13) 


n 

These  three  equations  for  the  Fresnel  zones  correspond  to  Eqs.  (B-l), 
(B-2),  and  (B-3)  for  the  irradiated  area. 

Consider  a numerical  calculation  of  the  dimensions  of  the 
insonified  area  and  the  dimensions  of  the  first  three  Fresnel  zones. 


The 


source  is  placed  6Q.Q  in.  (R^y)  from  the  surface  and  the  receiver  is  main- 
tained in  the  specular  direction  at  the  distance  of  6Q.0  in.  ( K,,^ ) from 
the  surface.  The  beamwidth  of  the  projector  (source)  is  9 deg.  The  results 
are  depicted  in  Fig.  &J*  for  a frequency  of  $5,6  kHz  for  grazing  angles  20 
and  bo  deg.  The  calculations  show  that  the  dimensions  of  th-  insonified 
area  are  slightly  larger  than  the  dimensions  of  the  first  Fresnel  zona. 

The  scattering  surface  boundary  shows  the  relative  size  of  the  model  surfaces 
to  the  ineenified  area. 

If  the  be-uswldth  of  the  source  la  varied  from  6 deg  to  Id  tfcg  for 
the  situation  just  described,  the  inson;  fled  arc a will  enclose  slightly 
less  than  one  Fresnel  sons*  to  almost  tv<  Fresnel  zone*. 

It  is  difficult  to  ascertain  with  analytical  precision  ih*  part 
of  the  scattering  surface  that  should  b*  taken  as  the  active  scattering 
region.  However.  sine«  "nly  a few  Fresnel  rones  are  tnsonified  in  the  expert* 
mental  study  in  this  paper,  the  insonifird  area  will  b**  taken  ns  the  active 
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FIGURE  B*4 

THE  FIRST  THREE  FRESNEL  ZONES  AND  THE  INSONIFIED  AREA 
FOR  A BEAMiftDTH  OF  9 d«9 
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scattering  region.  Thus,  the  limits  of  the  integrals  appearing  in 
Chapter  II  and  III  will  he  taken  as  the  values  of  A^,  and  l calculated 
from  Eqs.  (B-l),  (B-2),  and  (B-3),  respectively. 

It  is  noted  that  the  approximate  dimensions  of  the  insonified 
area  A^,  A^,  and  I are  used  as  the  limits  rather  than  the  exact  dimensions 
of  the  insonified  areas  a and  b.  This  approximation  is  equivalent  to  inte- 
grating over  a rectangular  area,  which  approximates  closely  the  elliptical 


Table  B-I  lists  the  values  for  A^,  A^,  and  l for  grazing  angles 
from  6 deg  to  80  deg  for  the  previous  cited  numerical  example  for  a 9 deg 
beamwidth  of  source.  The  dimensions  A^,  A^,  and  t are  calculated  from 
Eqs.  (B-l),  (B-2),  and  (B-3),  respectively.  The  irradiated  surface  area 
lor  grazing  angles  from  6 deg  through  20  deg  is  calculated  from  the  relation 

2i(A1  + A2) 

whereas  the  area  for  grazing  angles  from  22  deg  degrees  to  Qo  deg  is 
calculated  from  Eq.  (B-6). 
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DIMENSIONS  OF  INSONIFIED  AREA  FOR  VARIOUS  GRAZING  ANGLES 


Grazing 

Angle 

^i 

^2 

1 

Insonified 

Area 

Degrees 

Inches 

Inches 

Inches 

Sq.  Inches 

6 

16.00 

16.00 

4.72 

302.21 

8 

16.00 

16.00 

4.72 

302.21 

10 

16.00 

16.00 

4.72 

302.21 

12 

16.00 

16.00 

4.72 

302.21 

14 

14.84 

16.00 

4.72 

. 291.22 

16 

13-44 

16.00 

4.72 

278,06 

18 

12.30 

16.00 

4.72 

267.28 

20 

11.35 

16.00 

4.72 

258.32 

22 

10.55 

15.65 

4.72 

211.05 

24 

9.07 

14.10 

4.72 

190.06 

?6 

9.28 

12.84 

4.72 

175.55 

28 

8.76 

11.81 

4.72 

159.72 

30 

8.31 

10.93 

4.72 

148.38 

52 

7.91 

10.20 

4.72 

138.30 

34 

7.56 

9.56 

4.72 

150.61 

36 

7.25 

9.01 

4. 72 

123.53 

3« 

6.97 

8.53 

4.72 

117.36 

40 

6.72 

8.11 

4.72 

111.94 

42 

6.49 

7.73 

4 ,73 

107.15 

44 

6.?9 

7.40 

4.7? 

102.8a 

Li> 

6.10 

7.10 

4.72 

99.10 

5 .v|* 

6,84 

4.7? 

95.71 

50 

6.60 

4.72 

9f-.66 

J i 

$.6$ 

6.39 

4.72 

3 4 

5.5? 

649 

4.71 

#7.4$ 

''*> 

$.41 

6.0? 

4.7? 
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£0 
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#1.39 
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4.71 
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$.S6 
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4.94 
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4.71 
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nm 
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nM 
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TABLE  B-X 

DIMENSIONS  OF  XNSONX5TED  AREA  FOR  VARIOUS  GRAZING  ANGLES 


Grazing 

Angle 

*1 

*2 

I 

Insonified 

Area 

Degrees 

Inches 

Inches 

Inches 

Sq.  Inches 
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6 

16.00 

16.00 

4.72 

302.21 

8 

16.00 

16.00 

4.72 

302.21 

10 

16.00 

16.00 

4.72 

302.21 

12 

16.00 

16.00 

4.72 

302.21 

14 

14.84 

16.00 

4.72 

291.22 

16 

13.44 

16.00 

4.72 

278. 06 

16 

12.30 

16.00 

4.72 

267.28 

20 

11.35 

16.00 

4.72 

258.32 

22 

10.55 

15.65 

4.72 

211.05 

24 

9.87 

14.10 

4.72 

190.06 

26 

9.28 

12.84 

4.72 

173.35 

28 

8.?6 

11.81 

4.7? 

139.72 

50 

8.31 

10.93 

4.7? 

1.48.38 

52 

7.91 

10.20 

4.72 

138.80 

54 

7.56 

9.56 

4.72 

150.61 

5° 

7.25 

9.01 

4.78 

123.3? 

3« 

6.9/ 

8.53 

4.7? 

117.36 

uo 

fc ./? 

6.11 

4.72 

III.94 

42 

6.49 

7.73 

4.72 

107.13 

44 

6.?9 

7.40 

4.7? 

I0f.se 

4 ,6 

6.10 

7.10 

4.72 

99.10 

48 

5 .*'5 

6.84 

4.7? 

95.71 

;o 

!.?■* 

6.60 

4,78 

gf.il 

-,s 

5.65 

6.39 

4.78 

8*90*2 

5’4 

§.5? 

6.19 

4.7S 

S7,4§ 

Y-> 

5.41 

6.08 

4 71 

I5.8S 

5.51 

5.66 

4.7? 

33.21 

to 

3.22 

3.71 

U.7g 

-31.  $* 

545 

3.56 

4.7® 

79.76 

<>4 

5.06 

5.46 

4 .72 

1^.86 

fei> 

4.<*i 

5.36 

4.78 

76.# 

v# 

4.94 

5.86 

4.?,? 

”5 . /t 

V 

4.40 

5.17 

4.78 

74.JJ 

.*;■ 

4.«4 

5.10 

*.78 

■ ?*.*? 

i'4 

-.50 

3.0$ 

4.72 

72.# 

4.77 

4.p 

4.7? 

7S.8& 

•a 

4.7J 

4.«a 

6.7® 

; 71 .6® 

4./J 

4.86 

4.T2 
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EVALUATION  OF  THE  PRESSURE  REFLECTED  FROM  A PLANE  SURFACE 


In  Chapter  II,  Eq.  (17),  we  encountered  the  integral 


ik  7 r i [zi  z?  ik(Ri+Rs) 

Po  - J J R-R-  lR-  r R*  e 

A 


dx.jdy 


(C-l) 


where  the  limits  A-,,  and  I are  the  dimensions  of  the  acti\e  scattering 
region  defined  in  Appendix  B.  Calculations  shall  be  limited  to  the 
specular  direction. 

The  Fraunhofer  approximation  is  the  neglect  of  the  quadratic’  : 
and  higher  order  terms  in  the  expansion  of  the  term  R^  + R’  appearing  i-n  •' 
the  expedient ial  of  Eq.  (0-1).  For  the  specular  direction  and  for  the 
Fraunhofer  approximation  Eq.  (A-l4)  in  Appendix  A reduces  to 

Ri  + * R10  * R?0 


Setting  z.  « IL-  sin  i<  , ? - R„'*in  4 and  H!  o R 

A J V « J £V  i 

“■£  * R;>o  ln  denominator  of  Eq.  (C-l),  one  obtains 


l "10 


* t { A,  .♦  i\, ! sin  i iji  ♦K(S,Q+R^)j 


(c-») 


v ■'  ^ 
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•where 


* « . v:  i H tf 


Iv  §V 


% * i * Wa  * 


0 eO 


% la  the  asplUv-a*,  and 
Is  the  phase. 
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(C-$) 


( C-*» ) 


.■fiiiWSgL 


1 
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M 

Mi 


The  pressure  squared  P0P0>  which  is  proportional  to  the  intensity  reflected 
from  a plane  surface,  is 


op  p p 

* + ^}-  sin 


p p = 
*0  0 


2 2 2 
n R^q  R20 


(C-5) 


For  the  Fresnel  approximation  the  quadratic  terms  are  retained 
in  the  expansion  of  + R^.  The  desired  expansion  given  by  Eq.  (A- 19) 
in  Appendix  A 

«i  * "a  a *io  * ho  ♦ * <xi  sl“2  * * & ■ 

reduces  Eq.  (C-l)  to 


, \ |x?«tn8*l  | 

, IK  .u.  * ‘k<  W 7 ,upx~L  r 1 1 
*.'SSg£'  J 6 M6  * 

-I 


. (C-6) 


Tie  integrals  can  not  be  integrated  in  closed  fora  but  yield  infinite 
as  ■•two.  They  are  culled  Fremiti  integrals  and  are  defined  as 


/ 003  (4  «S)du 


f r 9 

3(u)  «•  ; sis  u )du 
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The  amplitude  is 


hw : %)  • sa  * %ag  * Ins  4 %)  * m\  : %i]gT  1 

(R10  + R20‘> 


and  the  phase  is 


CA  + C„ 

h 

sA  + sA  r 

, ^3 

Ls< 

C.  + C.  1 

[h  *t 

+ c 

l 

S.  + S. 

All 

<&Q  = « + ^{\Q  + R2q)  - arctan 


The  pressure  squared  is 

'iv%i-**iv%r*  Nv%rc*tvsi 


p p W -» 

0-0 


3 


(R10  - R2o)c 


(0-10) 


(u-11) 


(c-12) 


UH 


i 


rm 


The  amplitude  and  phase  obtained  in  the  Fraunhofer  approximation 
can  be  obtained  fro©  the  amplitude  and  phase  derived  in  the  Fresnel 
approximation  if  it  is  assumed  that  the  argument  of  the  Fresnel  integrals 
is  very  small  so  that 

C(u)  58  u , and 

• S(u)  « 0 . 

Physically,  this  approximation  implies  that  both  tiie  source  and  receiver 
oust  be  moved  to  infinity  with  the  restriction  that  the  source  have  a very 
narrow  beam  spread  so  that  the  dimensions  of  the  insoaified  area,  , .\s, 
and  ! are  very  small  compared  to  the  distance  from  the  surface  to  the 
source. 

these  statements  assert,  that  it  is  theoretically  impossible  to 
predict  the  intensity  reflected  from  a plane  surface  In  the  specular 
direction  in  the  Fraunhofer  approximation,  h comparison  of  sqs.  (C-b) 
and  (C-ll)  will  substantiate  the  invalidity  of  the  fraunhofer  approximation 


’ -s 


JyyA,,: 
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in  predicting  the  correct  phase.  One  would  expect  ■'  phase  change  at 
the  receiver  relative  to  the  source  of  « radians  at  the  surface  and 
k(\o  + R2cP  radians  as  a result  of  the  travel  path.  Figure  C-l  depicts 
the  phase  in  the  Fraunhofer  and  Fresnel  approximation  for  the  case 
where 

f =95.8  kHz, 
k « 10.323  in."1, 

R^O  = R20  = iu* 

The  values  for  A^>  Ag,  and  i for  a 9 deg  beamwidth  of  source  are  given  in 
Appendix  3.  The  term  « + k(^Q  + 11 ' not  inc^ude<i  iu  tiie  calculations. 

For  this  numerical  example  the  Fresnel  plume  lies  between  20  and  30  deg 
for  grazing  angles  greater  thar.  30  deg.  Tin?  Fraunhofer  phase  is  -90  deg 
for  all  grazing  angles. 

The  amplitude  in  the-  Fraunhofer  approximation,  Eq.  (C-3)»  and 
the  amplitude  in  the  Fresnel  approximation,  Eq.  (C-10),  are  plotted  for 
this  numerical  example  on  a relative  “scale  for  various  grazing  angles  in 
Fig.  C-2.  Over  a wide  range  of  grazing  tingles  (kO  to  80  deg)  the 

Fraunhofer  amplitude  is  approximate!:/  t.i  tinea  the  Fresnel  amplitude. 

* 

The  term  p p that  is  proportional  to  the  intensity  in  the 

J Q £ 

Fraunhofer  approximation,  Fq.  (C-5),  and  the  tew  pp  in  the  Fresnel 
approx  heat ion,  gq.  (C-l«)  are  platted  for  this  numerical  example  in 
Fi|.C«5.  Tte  Fraunhofer  intensity  is  aver  Id  times  the  intensity  in 
the  Fresnel  approximation  for  grating  angles  greater  than.  §0  deg. 
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FIGURE  C-l 

RELATIVE  PHASE  IN  THE  SPECULAR  DIRECTION 
PLANE  PRESSURE  RELEASE  SURFACE 
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INTEGRATION  OF  TERM 


In  Chapter  III,  Eq.  (68),  we  encountered  the  integral 


v1*  ZtA  • 11  , ^2  1 A2  £ if/  tv2  , ,,2 

k (Fq)  sin  y p p p r --  (x1-x1)  +( y-y  ) 

Jl " 2 r2a2  J J J J e 

10  20  o -I  -£ 


2 _ 2 .2 


X cos 


fk  . 2 ,1  2 ,2\  k . 2 ,2 

Li  s“  *fi  - xi.  + I (y  - » 


)j  d^dydxj 


dydx.'dy' 


Making  the  change  in  variables  defined  by  Eq.  (58),  one  finas  for 


A /_2v  . b A i A i iC  f>2 

k (fo)  sin  v r r r r (u_u  } ' 

_ —"5  2 ~ 2'  / / / e 

2*"TU  R_  A"  0 - J 


r+Cy-y1)2 


10  20  c -A  -£  -A  -£ 


x cos  | - u'Ain2  V + | (-'V  - ^)(u  - u^sin*"  0 

+ | (y~  - y*4*)  dudu’dydy* 

Introducing  relative  and  center  of  gravity  coordinate.-  defined  by 
gqs.  (61)  and  (62),  one  obtains 


, tG"- 


. , , A £ 2(A-X  ) 1/ 

h 6;)  ntu*  v rr  r ° f 

<S  * — !?-— ?r~“  / 1 ' Is 

* H ' 8 * A & v v w , » . 

tu  20  0 e o •,?(  A-x^)  -S(l-y^) 


- - »\  ^ “ 
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l\y 
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2k  2 A A 
■sr-  y dx  dy  ' 
R ',oj  o 


dy^dfcdn 


l\t  + 


A. 

_i  p + — 

2 i 6 R 


§ — p cos  <p 
n a p pin  <p 


-p/a 

e ecu? 


M * |“  (sir.  y) 


pdad<? 


A A, 


nev,  eet 


%'q  efetsvin 


U a TT*  Jf 

S e- 


H a P eea  y 
ih?  g in  y 


«©/ Si 

<3  ' 1 «*©S 


*/ 

e p 
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» / H'  ♦ JS  <s«8(?P  - >J  ;'4’4” 
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The  term  (cp  - y)  ranges  from  - 7 to  2n  - 7 as  cp  goes  through  its  range 
of  o to  2n;  the  integrand  assumes  all  the  values  that  it  would  if  the  value 
of  7 were  zero.  Thus, 


CJt  “ 

= f f e p^a  cos  p J + N 


o o 
o o 


cos  cp 


pdpdcp 


From  Gradshteyn  and  Ryzhil.  (integral  3.715>-19>  P*  ^*02)  one  finds  for  the 
integral 


f r ^ 2 2 f~2 

i cos  p J M + N cos  #\  dcp  - 2nJo  pj  M + N"  J , 


where  J is  the  2eroth  ordered  Bessel  function.  The  expression  for  I 
o 

becomes 


I * 2n 


f e-o/"  oj 

/ 2 2 

p J M + K 

c 

o 

V 

d*3 


Using  an  integral  solution  (6.623-2,  p.  712)  hi  Gradshteyn  and  Ryshik, 
one  obtains  for  I 


I ““ 

a 


it  sin  ; 1 

f \ 

V + __  _ 

v * J 

1 ° * « 

.-hi  eh  gives  for  J 
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; K rir/  tVI  / ^ 

“■4 1 — r — 
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Integrating  over  y q by  employing  the  table  of  integrals  by  Gradshteyn 
and  Ryzhik  (integral  2.271-5,  p.  86),  one  gets 


■ . j • . •;  : - A rj  Jgfwrt ■ . •.;••  • Vft;  t.V- 1 - V.  Ij-'v. 

tSUafcVSL  ■’■* :l  ■■JliiW-w owintW^ 


Integrating  over  yQ  by  employing  the  table  of  integrals  by  Gradshteyn 
and  Ryzhik  (integral  2.271-5,  p.  86),  one  gets 


, ..  *lxdi 

0 ~ p 2 2 2 

A sin  y 

10  dv  0 
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to  get  for  J+ 
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Prom  this  equation  one  can  write 


u (Fo>  f dv 
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and  since 


one  obtains 
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The  substitution  of  the  expressions  for  A and  B,  Eq.  (D-9)>  &nd  the 

expression  for  A , Eq.  (C-10)  of  Appendix  C,  into  J.  gives 
o «•» 
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* <0  •-* 


J.  3 


r.  'l' 
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2/  -1 


,(here  f and  C are  the  Fresn*'1  integrals  defined  by  Eq.  (C-9)  Appendix 


For  the  special  case  where 


S and  C both  approach  one-half,  and  Eq.  (D-ll)  reduces  to 
2 / 2\  . 2 
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4k  sin  \{i 


X arccot 


2ka£A  sin  ^ 
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2 2 4 

i " + A sin  xir  + 
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53? 
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which  can  he  compared  with  the  terra  obtained  by  Gulin  (1962),  the 
result  of  using  a Gaussian  covariance  function. 
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( D— 12) 


DISCUSSIOH  OP  J2  TERM 


In  Chapter  III,  Eq.  (69),  we  encountered  the  Integral 


(E-l) 


Making  the  change  of  integration  variables  as  depicted  in  Fig.  t),  the 
author  could  not  find  an  analytical  expression  for  j,,,  but  approximate 

C 

expressions  showed  to  be  very  small  compared  to  . 

The  conclu,  >,on  that  is  much  less  than  .7^  is  borne  out  by 
examining  the  work  of  Gulin  (19o§)  who  has  evaluated  the  integrals  and 
Jg  for  the  Gaussian  form  of  the  covariance  function,  liuaerieal  evaluations 
of  Ms  ‘,xpr«s8io89  are  illustrated  in  Fig.  g-1  for  the  case 

RJ0  ' rV<i  * ln’ 


«.  Ov\*)l  in. 
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k a 10.  J5  in‘- 
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C.SQd  insh.  ths  «5*d  far  a&i  t *fr  listed  is  fafele  8*1 

A * 

in  Appendix  S.  At  a feyleiPh  p&rsmbt**r  at  0 J (tfee  theory  t«  valid  only 
fey  Enylelgh  parameters  less  thaws  one'  the  xaplitusi?  fluctuations. 

^ * Jg)  ■ *\  exesed  the  pi'^tse  flust uatissSs.  - Ug,}  * by  apftf'bxijs&fely 
b percent.  §lftr©  for  small  Saylei^h  parsaet^rs  the  ^plitude  and  phase 
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FIGURE  E-l 

THEORETICAL  FLUCTUATIONS 
IN  THE  SPECULAR  DIRECTION 
FOR  A GAUSSIAN  COVARIANCE  FUNCTION 
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fluctuations  are  approximately  of  the  same  magnitude,  Gulins'  expressions 
for  and  show,  for  the  example  cited,  that 

J2  « Jx 

Figure  E--2  shows  that  the  term  for  the  case  of  an  exponential 
covariance  function  does  not  differ  much  from  the  term  for  the  case 
of  an  Gaussian  covariance  function,  especially  in  the  region  of  Rayleigh 
parameters  less  than  one  where  the  theory  is  valid.  Thus,  one  concludes 
that  the  Gulin  result  of  Jg  « for  the  Gaussian  case  is  applicable 
also  to  the  exponential  case. 
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